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ABSTRACT 


A high temperature stage fitted to a Siemens Elmiskop 1 electron microscope 
has been used to study directly the annealing of prismatic dislocation loops 
in thin foils prepared from quenched and aged bulk specimens of aluminium. 
In the temperature range 170—200°c the loops are observed to shrink and 
eventually to disappear. Observations on dislocations other than loops 
also present in the thin foils show that they become smoother while the loops 
are disappearing. The results are interpreted in terms of climb of dislocations 
by vacancy emission, and in the case of loops in thin foils, where vacancy 
supersaturation effects are thought to be negligible, the climb rate is governed 
by the line tension and curvature of the loops. The radius 7 of a loop as a 
function of time ¢ is given by the theoretical expression r=r,[ 1 —(¢/7)]*/2 where 
7 is the lifetime of a loop of initial radius 7). The temperature dependence 
is primarily t=7, exp (H/kT) where #H is the activation energy of self 
diffusion. Comparison with experiment gives estimates of the activation 
energy in the range 1-2—1-3 ev. 

Annealing experiments with bulk material have been carried out in parallel 
with the above experiments. They show that the loops anneal out in the 
same temperature range. The details are different however owing to vacancy 
supersaturation effects. Larger loops appear to grow at the expense of smaller 
ones. The theory gives a qualitative explanation of the annealing of 
electrical resistivity in quenched aluminium observed by other workers, 
although the shape of the isothermal annealing curve is difficult to predict 
quantitatively, owing to its dependence on the complex annealing behaviour 
of loops in bulk material. The experiments also suggest that the temperature 
rise of thin foils examined in the electron microscope due to heating by the 


electron beam is very small. 


§ 1. INTRODUCTION 
WHEN a metal is rapidly cooled from high temperatures large numbers of 
vacant lattice sites in excess of the equilibrium concentration may be 
quenched-in. The changes in physical properties produced by quenching 
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and the variation of these changes on subsequent annealing are therefore 
interpreted in terms of the properties and behaviour of vacancies (Maddin 
and Cottrell 1955, Bauerle and Koehler 1957, Koehler et al. 1957, Panseri 
and Federighi 1958, Federighi 1959, DeSorbo and Turnbull 1959, Kimura 
et al. 1959a,b, Meshii and Kauffman, 1959, for review see Cottrell 1957). 
An investigation of the annealing out of the increase in electrical resistivity 
of aluminium produced by quenching has recently been carried out by 
Panseri and Federighi (1958) and by Federighi (1959). These workers 
showed the presence of two stages in the isochronal annealing curves of 
specimens quenched from temperatures greater than ~400°c. ‘The first 
stage took place in a few minutes at room temperature (Wintenberger 
1957, Bradshaw and Pearson 1957, DeSorbo and Turnbull 1959) and was 
associated by Panseri and Federighi with the movement of vacancies to 
form clusters of some unspecified type. The second stage occurred at 
temperatures at which self-diffusion is important and was therefore 
associated with the dispersal of the clusters. 

Recent observations (Hirsch, Sileox, Smaliman and Westmacott 1958) 
using the transmission electron microscope technique, have revealed the 
presence of small prismatic dislocation loops in thin foils of aluminium 
prepared from quenched and aged bulk material. These are thought to be 
produced by the condensation of excess vacancies into disc-like cavities 
which subsequently collapse to form loops of whole dislocation line. How- 
ever, aS was pointed out by Federighi (1959), the identification of these 
loops with the clusters postulated by Panseri and Federighi was not 
unambiguous. Thus it was possible that while a fraction of the vacancies 
formed dislocation loops, another fraction formed clusters so small that 
they are not detected on transmission electron micrographs. Stage II of 
the resistivity annealing curves could then be associated with the annealing 
of small vacancy clusters and not necessarily with the disappearance 
of loops. It is therefore necessary to decide this point by experiment. 
Two types of experiment may be carried out using the transmission electron 
microscope. In, the first a high temperature stage in the microscope is 
used to observe the annealing out phenomena directly in the thin foils. 
The behaviour of the loops on annealing can therefore be studied in, detail. 
In the second, the annealing treatment is carried out on bulk material 
which is subsequently thinned and examined in the electron microscope. 
This paper describes the results of both types of experiment for aluminium 
and discusses their physical significance. 


§ 2. EXPERIMENTAL TECHNIQUE 

Aluminium of Johnson-Matthey ‘specpure’ quality, nominally 99-99% 
pure, and in the form of foil 0-003in. thick, was quenched from about 
600°c by dropping under gravity from a vertical tube furnace into iced 
brine. The foil was then allowed to age at room temperature for several 
hours. ‘The foils were thinned by electropolishing in an 80% ethyl alcohol 
—20% perchloric acid bath (Tomlinson 1958). Observations were 
recorded on photographic plates in the usual way, and also in some of the 
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high-temperature stage experiments by direct ciné-photography of the 
fluorescent screen of the electron microscope using equipment described 
by Whelan e¢ al. (1957). For the bulk material experiments a quenched 
and aged foil was cut into two parts. One part was thinned and examined 
directly to establish the existence of loops. The other part was annealed in 
an oil bath for six minutes prior to thinning and examination. Annealing 
temperatures of 100, 160, 170, 180, 190, 200, 210 and 220°o were employed. 

Examination was carried out with a Siemens Elmiskop 1 electron 
microscope operating at 100 kv. 

Figures 1 (a), (6) and (c) show the main, details of the high temperature 
stage used for annealing the thin foils during observation. The specimen 
holder (fig. 1 (a) and (b)) is similar in shape to the conventional holder and 
is inserted through the airlock into the bore of the objective pole-piece in 
the usual manner. The specimen holder contains a two-pin socket E, 
insulated from the body of the holder and carrying a filament type speci- 
men support AB. This support consists of two 0-4 mm diameter platinum 
legs B spot-welded to a 200 mesh stainless steel grid A and secured in 
the socket by two screws D. The grids are cut from RCA } in. diameter 
stainless steel specimen supports. The spacing between the platinum 
legs at the grid is 1-5 mm, and the width of the grid is such that seven 
strands of wire conduct the heating current. Direct current is passed 
through the grid by means of a two-pin plug A at the end of a long rod 
inserted through the hole in the objective lens normally occupied by the 
stereo-drive (fig. 1 (c)). With this equipment therefore it is not possible 
to tilt the object during observation. The rod can move to follow the 
stage movement and can be withdrawn to an inoperative position when 
the specimen is removed. The instrument also requires a modified stage 
carriage (not shown) similar to the conventional carriage but without the 
stereo-mechanism. More complete details of this equipment have been 
reported by Whelan (1958). 

The specimens are mounted by removing the end cap C in fig. 1 (a). 

. The grid is dipped in a dilute solution of polybutene in xylol, which acts 
as an adhesive ; excess solution is blotted off and the specimen is placed 
in the centre of the grid as in fig. 1 (0). 

The compact design of the stage prevents measurement of the tempera- 
ture of the grid during observation. However grid temperature versus 
heating current curves have been determined by calibration of individual 
specimen supports in a high vacuum system. A Pt, Pt-13% Rh thermo- 
couple of 48 s.w.g. (supplied by Johnson-Matthey) was placed as closely 
as possible to the centre of the grid. Millivolt readings were converted 
to degrees centigrade using Johnson-Matthey conversion tables. Measure- 
ments were reproducible to about +2°c at 170°c. The current required 
to heat the centre of an average grid to 200°c was about 1-1 amp, the 
working temperature being attained a few seconds after switching on. 
The variation of temperature over the grid was not investigated. How- 
ever, in a few cases a grid was recalibrated with a new thermocouple 
after use. The temperature usually differed slightly from the previous 
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High temperature stage for the Elmiskop 1. 
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estimate. From such measurements the inaccuracy in stage temperature 
obtained from heating current in this way is very roughly estimated to 
be about +10°c. This does not include of course any rise of temperature 
of the thin foil due to heating by the electron beam. To minimize this 
last effect the double condenser was used with an illuminating spot about 
10 » diameter, with a 600 . diameter aperture in condenser 2, and with a 
total beam current ~ 15 pa. 


§ 3. EXPERIMENTAL RESULTS 
3.1. High Temperature Stage Experiments 


When a specimen containing loops is heated in the electron microscope 
the loops are observed to shrink and eventually to disappear completely. 
The loops begin to shrink at a noticeable rate at a stage temperature of 
about 170°c; the rate of shrinking increases rapidly with increasing 
temperature. For example, timed sequences of micrographs and ciné- 
sequences have shown that the rates of decrease of radius of a loop about 
250A radius are about 0-2, 1-5, and 33 Asec—! at temperatures ~ 170, 200 
and 220°c respectively. This shrinking is quite distinct from the process 
of prismatic slip observed by Hirsch, Silcox, Smallman and Westmacott 
(1958). Figure 2 (a), (b), (c) and (d) (Pl. 1) is a sequence of micrographs 
taken, on photographic plates at an, estimated stage temperature of 187°C. 
The time interval between successive micrographs was not recorded but 
is estimated to be about one minute. Most of the loops in fig. 2 (a) 
diminish in size and disappear during the sequence ; loop A remains through 
the whole sequence although diminished in size. In general in the thin 
foils, loops which are originally round remain round as they shrink, i.e. 
there is no tendency to form angular shapes. Also all the loops (large and 
small) appear to shrink in the thin foils. There is no tendency for large 
loops to grow at the expense of small ones. 

The rate at which a loop shrinks is not uniform but is a function of its 
radius. A small loop shrinks much faster than a large loop, i.e. the rate 
of shrinking of a loop is accelerated as it grows smaller. Figure 3 (PI. 2) 
is a sequence of micrographs taken on photographic plates of a loop 
shrinking at a stage temperature of 197°c; the times in seconds at which 
these exposures were made are marked on each micrograph. The loop lies 
on an inclined plane, and the distance between the more steeply inclined 
and less visible parts of the loop (Hirsch, Howie and Whelan 1958, 
1959) has been taken as a measure of the loop diameter. Figure 4 shows a 
plot of loop radius versus time. There is some uncertainty in the exact 
position of the centre of the dislocation line as judged from the contrast 
and this is indicated by the vertical spread in the points on fig. 4. How- 
ever, the relative positions of the points are probably more accurate. With 
single exposures on photographic plates it is not possible to record the 
shrinking at very small radius owing to its rapidity. However using 
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ciné-techniques the steep slope of the radius-time curve at small radius has | 
been verified, although it is difficult to determine its functional form 
with accuracy. 

The observations suggest that the loops are acting as vacancy sources, 
i.e. that the vacancies which condensed to form the loops after quenching 
are dispersing. Assuming that the only force acting to aid climb is that 
due to the line tension and curvature of the loop, the theory given in detail 
in §4 suggests that the radius time curve of a shrinking loop should be 
approximately a parabola. Accordingly such a curve has been fitted to 
the experimental points of fig. 4 and the points are seen to be consistent 
with this curve. The dislocation loops are therefore climbing by vacancy 
diffusion, in a sense to annihilate themselves. These results may be taken 
as direct proof of the hypothesis of climb of dislocations. Figure 5 (a), (0), 
(c) and (d) (Pl. 3) shows a sequence of micrographs from a region 
containing both loops and irregular dislocations (e.g. fig. 5(a) at A). 
During the sequence, while the loops are disappearing, the irregular 
dislocations assume a noticeably smoother shape (e.g. at A). This effect 
has been observed a number of times and is thought to be due at least 
partly to climb. 


Radius versus time plot for the loop in fig. 3. A parabola has been fitted to 
the experimental points. 


Relatively little slip of the loops has been observed. Occasionally how- 
ever a loop has been observed to jump suddenly as a whole over a small 
distance, presumably along its slip cylinder, after it has shrunk to a very 
small size before it finally disappears. The explanation of this effect is: 
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not known with certainty ; however prismatic slip is a very rapid way of 
transferring the last remaining vacancies as a whole to the surface as a 
nearby sink. 


3.2. Bulk Material Experiments 


The density of loopsin the quenched and aged specimensis ~ 6 x 1014em-3 
assuming a foil thickness of ~ 10004. No noticeable change was observed. 
in thin foils prepared from bulk material which had been annealed for six 
minutes at temperatures below 180°c. In specimens annealed at 180°C 
the density decreased slightly to ~1:3x10!em-*. A reduction in loop 
density was quite marked in foils annealed at 190°c and 200°c in which the 
densities were about 5 x 1018 and 5 x 101°cm- respectively. These figures 
should be regarded as order of magnitude estimates only since variations 
were observed from specimen to specimen, and from point to point on any 
one specimen. For reliable estimates of loop densities the foil thickness 
must be determined (for example from the projected width of a slip trace 
left by a dislocation moving on a {111} plane). It has not always been 
possible to determine the foil thickness in this way, and an average thick- 
ness ~1000A has been assumed. Experience has shown that the foil- 
thickness in areas transparent to electrons is usually of this order. The 
change in number of loops on micrographs from as-quenched material and 
from material annealed at 190°c and 200°c is too large to be attributed 
to differences of foil thickness. 

Figure 6 (a) (Pl. 4) shows a micrograph of a foil prior to annealing 
and fig. 6(b) (Pl. 4) and fig. 6(c) (Pl. 5) show micrographs of foils 
annealed at 190°c. Several large loops can be seen in these micrographs. 
This increase in size accompanying the decrease in loop density is also 
apparent in the micrograph of a foil annealed at 200°c shown in fig. 6 (d) 
(Pl. 5). In specimens annealed at 210°c and higher temperatures no 
loops were observed. The conclusion drawn from these results is that 
in bulk material large loops grow at the expense of smaller loops; even- 
tually however they too either shrink and disappear or become large 
enough to be incorporated in the dislocation network of the crystal. 

Sometimes the larger loops left in bulk material after annealing have 
angular shapes (e.g. figs. 6() and (c) at A). Such loops may be parallelo- 
grams with edges in close packed directions although no experimental 
confirmation of this has been attempted. The shape of a large loop which 
has grown from a small one by vacancy absorption presumably depends on 
the detailed mechanism of growth of a loop (see § 4.1). 


§ 4, THEORETICAL INTERPRETATION 
4.1. Olimb of Dislocation Loops 


The observations leave little doubt that climb of the dislocation loops is 
occurring and it is therefore convenient to discuss the various factors which 
govern the climb of a prismatic loop. Friedel (1956) has considered the 


8 J. Sileox and M. J. Whelan on the 


climb of a dislocation due to vacancy diffusion in terms of the movement of 
jogs. He gives the following expression for the velocity of a jog along the 
dislocatiou line: 


w nana Oo fee Hie Ea Pob\ _ Lay ; 
oe ara" exp | kT }L exe Po) \ () 


where z~ 11 (of the order of the atomic coordination number of the lattice), 
va is an atomic frequency (of order 1018); is the angle between the Burgers 
vector b and the dislocation line; E; is the formation energy and Hm the 
activation energy of movement of a vacancy ; H = #;+ Hm is the activation 
energy of self-diffusion; k7’ is the usual Boltzmann, temperature factor ; 
while F, and Fs are forces per unit length acting on the dislocation line 
tending to make it climb out of its slip plane. /, is a force due to an 
applied stress or due to the stress fields of other segments of dislocation. 
F', is a chemical stress due to a deviation of the concentration of vacancies 
from the equilibrium value at temperature T. /’s is equal to 

kT 

“ (2) 
where c and ¢, are the actual and equilibrium concentrations of vacancies 
at temperature 7’. 


Life 


log ES 


Fig. 7 
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Climb of an edge dislocation by jog movement (following Friedel 1956). 


Equation (1) may be given a simple physical interpretation. Consider 
an edge dislocation containing a jog in a simple cubic lattice as in fig. 7; the 
extra half plane of atoms is shown shaded. Under the action of a suitable 
force F', tending to make the dislocation climb the jog will emit vacancies 
and will move to the right with an average velocity given by the left-hand 
term in the square brackets of eqn. (1) together with the factor outside the 
brackets. However vacancies will tend to recombine with the jog which 
will therefore move to the left, and this is taken into account by the right- 
hand term in the square brackets. The resultant velocity of the jog is the 
difference between these two terms.. (The formula assumes that the 
force F’, does not affect the probability that a vacancy will combine with 
the jog, i.e. that the energy barrier to be surmounted when a vacancy 
combines with a jog is independent of the force F',.) Clearly the rate of 
climb of the dislocation normal to its slip plane in fig. 7 is 


V=¢ b/x 
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where x is the average distance between jogs. 6/2 is called the concen- 
tration of jogs c, by Friedel. ce is a dimensionless quantity equal to unity 
for an edge dislocation in a simple cubic lattice running at 45° to its 
glide plane. 

For a circular prismatic loop lying in a (111) plane in a face centred cubic 
metal such as aluminium the jog concentration will vary round the loop as 
shown schematically in fig. 8. ce will be of the order of unity at the 
segments of the loop lying between close packed direction in the (111) plane 
(i.e. at the corners of the inscribed hexagon in fig. 8). The jogs on these 


Fig. 8 


Circular dislocation loop in a (111) plane. The loop represents a disc of 
vacancies. The sides of the inscribed hexagon are parallel to close 
packed directions. The highest jog density occurs at the corners of 
the hexagon. Climb of the dislocation loops occurs by movement of 
jogs of opposite sign from adjacent corners and subsequent annihilation. 


segments will run along the loop in the directions indicated with arrows 
and will annihilate with jogs moving in opposite directions from the other 
joggy segments. Thus there is a tendency to produce a hexagonal loop 
lying inside the original circular loop as in fig. 8. However new jogs are 
expected to nucleate easily at the corners of the hexagon since this requires 
no additional line energy. The equilibrium form of the loop as it shrinks 
will be governed by the rate at which jogs are nucleated and the rate at 
which material is redistributed along the dislocation, i.e. on the rate of 
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diffusion between jogs. This is a theoretical problem involving jog kinetics 
analogous to the type of problem discussed by Frank (1958) and. Cabrera 
and Vermilyea (1958) for crystal growth kinetics. A similar theory could 
probably be applied to a shrinking loop, although the solution is not 
obvious. The experimental observations show however that as a circular 
loop shrinks it remains roughly circular, and this implies that jogs are 
nucleated and move at a rate sufficient to maintain a circular shape. It 
is worth noting here that in the converse process (i.e. when a small loop 
grows by absorption, of vacancies as in the bulk material experiments) the 
details may be different, as a little consideration of fig. 8 suggests. In 
this case new jogs will ultimately have to be nucleated by condensation of 
vacancies. This will involve an increase in dislocation, line energy and is 
therefore likely to be a difficult process. The dislocation loop will there- 
fore tend to have a low concentration of jogs, i.e. to be angular with sides 
parallel to close packed directions. This is probably the explanation of 
the parallelogram shaped loops visible in figs. 6(b) and (c). 

Detailed considerations of fig. 8 show that for a loop in a (111) plane of a 
face centred cubic lattice the maximum rate of change of radius of the 
loop is $ve, where ve is measured along close packed directions, i.e. Ce ~ $. 
In this theory for simplicity it will be assumed that the loops are circular 
and lie in planes perpendicular to the axis of the glide cylinder (i.e. in (110) 
planes). Thus all segments of a loop are edge dislocations (sinys= 1), and 
no dissociation into partial dislocations will occur. An average value of 
jog concentration cc ~ 4 will also be assumed to hold for all segments of a 
circular loop. Fora loop lying ina (111) plane the situation is not expected 
to differ much from this provided dissociation into partial dislocations can 
be neglected. The rate of shrinkage of the loop will then be given by 


et vee) ok ee 


where r is the radius of the loop at time f. 

It should also be mentioned here that eqns. (1) and (3) give the climb 
rate of a dislocation assuming that vacancy emission (or absorption) is 
the rate determining process. This will be the case if vacancies diffuse 
rapidly to sinks (such as free surfaces, grain boundaries, or other dislo- 
cations). If the diffusion is not sufficiently rapid a local supersaturation 
will be set up at the dislocation which will suppress the emission. The rate 
of climb in this case will be determined by transport of vacancies away 
from the dislocation by diffusion (the ‘saturation’ case discussed by 
Friedel 1956). The solution of this problem is complicated and depends 
on the geometrical conditions in any particular case (i.e. on the nature and 
distribution of sinks). This point is discussed in § 4.4 and in the Appendix. 


4.2. The Effect of Up-quenching 


During the annealing process the climb rate of a loop will be determined 
by the factor in square brackets in eqn. (1). It is assumed that both in 
the bulk material and in the thin foil experiments the metal is not under 
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any form of stress. In thin foils stresses may be produced by electron 
bombardment in the electron microscope; these however will be neglected 
since little prismatic slip was observed with low electron, beam intensity. 
Neglecting forces between loops, the only elastic force acting on a segment 
of loop will be that due to its line tension and curvature. This force tends 
to contract the loop. As far as the chemical stresses are concerned the 
experiments can be visualized as follows. The foil containing loops is 
assumed to have a concentration of vacancies equal to the equilibrium 
value at room temperature. The foil is then up-quenched from room 
temperature to about 180°c. This gives rise to an undersaturation of 
vacancies. Consequently the initial chemical stress is not equal to zero 
and acts in a direction tending to shrink the loops. A simple calculation 
shows however that to restore equilibrium only about one vacancy 
per loop is required (assuming a loop concentration of ~10!4cm-*) and 
equilibrium is rapidly attained, i.e. Fs=0. 

Once vacancies are being produced by the dislocation loops it becomes 
necessary to consider the rate at which they are absorbed by sinks. There 
will now be a supersaturation of vacancies and a chemical stress will build 
up in opposition to the line tension, stress, i.e. vacancies will tend to con- 
dense on the loops again. The discussion of this problem will be deferred 
to§4.4 At present it will be assumed that vacancies are removed by sinks 
faster than they are produced by loops so that the chemical stress is small, 
i.e. the right-hand term in the square bracket of eqn. (1) is very nearly 
unity. This is thought to be the case in the thin foil annealing experiments 
where the foil surface is a very effective vacancy sink (see $4.4 and 
Appendix). 


4.3. Climb of Loops in Thin Foils 


An approximate expression for the force F’, due to the line tension and 
curvature of a loop is (Frank 1950, Franz and Kréner 1955) 


ey 0 ioe {r/d) 
 4n(l—o) (rb) 
where the radius of the dislocation core is taken as b, Gis the shear modulus 
and a is Poisson’s ratio. (The core energy of the dislocation is neglected 
in this expression.) Putting fs=0 in eqn. (1) the radius r of a loop as a 
function of time ¢ is given by 


Pe0s 
a = —hzvab exp(— zm) | exp (Gr) 2 i Senge aes (1) 


With G=2:53x10" dynes cm? at 161°c, b=2:864, o=+ and 
k=0-86 x 10-4evdegreec™!, the reciprocal of the function in square 
brackets has been plotted as a function of r/b for temperatures of 161*c 
and 200°c in fig. 9. The curves approach infinity near r/b=1 but this 
corresponds to a range of r/b where the expression (4) is no longer valid. 
In the range of loop radius of interest therefore [exp Gch Li ig: 1] may be 
approximated by ab/r where «~72 at 161°C and ~ 64 at 200°c. dr/dt is 


F, (4) 
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therefore only slightly dependent on temperature through the factor « 
and depends mainly on temperature through the term exp (—/kT) 
in eqn. (5)7. 

Fig. 9 


r 4 
§(r) 


0) 50 100 150 200 250 


&(r) ——-—-— linear approximation 
Plot of the function 
ts Gb* log (r/b) -1 
fi) = [ ex (gEra=a (rib) ) -1] 


as a function of r/b. Note that the function is approximately linear 
over the range of 7/b of interest. 


ee e = —benubPaexp ( zn) = See 
Equation (6) can, be integrated to give 
roory(l tel easy) cee a 
where ry is the initial radius of the loop at t=0 
7r=7, exp (E/kT) wi) Ties Sete sane, Merl 
and where To= Tee] (zvab2x). VA sie S55 


+ Note added in proof.—Calculations show that an increase of the core radius 
of the dislocation has a similar effect on €(7) as an increase of temperature, 


e.g. at 161°c and 200°C changing the core radius to 2b causes « to decrease 
by about 10%. 
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Equation, (7) shows that the radius—time curve of a shrinking loop should 
be parabolic as shown in the curve fitted to the experimental points of 
fig. 4. These points were obtained from the loop in fig. 3 climbing at 
~197°c. In this case rp ~ 330 A, 7 ~ 138 sec, kT’ ~0-0405ev. Substituting 
these figures in eqns. (8) and (9) with z~11, va~ 104, a~ 64, H is found 
to be 1-29ev. The error in # due to 10°¢ error in temperature estimation 
is about 2%. The result also assumes the value of vg. The activation. 
energy depends logarithmically on this factor and is relatively insensitive 
to it. Errors due to rises in temperature of the foil in the electron beam 
are not known, but are thought to be small (see §6). Other estimates of # 
from the rate of shrinkage of certain, loops (for which complete radius—time 
curves were not obtained) gave values lying in the range 1-2 to 1-3ev. 
Although these values of the activation energy cannot be claimed to be 
accurate owing to assumed factors and to uncertainty in the accuracy of 
the temperature, the values are consistent with those obtained by other 
workers: Bradshaw and Pearson (1957) give 1-2ev; Federighi (1959) 
gives 1:33+0-05ev; DeSorbo and Turnbull (1959) give 1:3140-08ev; 
while Spokas (reported by Federighi 1959) has obtained 1-4 + 0-lev. 

The rates of decrease of radius of a loop of radius 250A can now be 
calculated from eqn. (6). Assuming a~ 70 for all temperatures of interest, 
and taking H = 1-32 ev (the average of the values quoted by Federighi and 
De Sorbo and Turnbull), the calculated rates are ~ 0-13, ~ 1:3, ~4Asec™} 
at temperatures of 170, 200 and 220°C respectively. Considering the errors 
involved in estimating temperature and in assumed factors these rates 
are in good agreement with the observed rates mentioned in § 3.1. 

The observations are therefore consistent with the theory of shrinking 
: of loops by climb. The temperature dependence of the rates is best 
illustrated by considering the lifetimes + obtained from eqn. (8) of a 
loop of radius ~ 1404. This is the average loop radius of the distribution 
shown in fig. 10, which is a histogram of loop sizes for a typical specimen 
quenched from 600°c (this histogram was obtained from fig. 1 of Hirsch, 
Silcox, Smallman and Westmacott (1958)). 7 has the following values 
for several temperatures (« has been taken ~ 120 at 20°c and ~ 70 for all 
other temperatures) : 


Bae 20 160 170 180 190 200 210 220 
7 230years 11-6min 52min 23min 1-:2min 35sec 17-5sec 8:8 sec 


It is seen that the loops will shrink at a noticeable rate for temperatures 
greater than ~170°c which is precisely the range in which the loops are 
observed to anneal out in the thin foils. 

Equation (6) can be written as follows (the plane of the loop is 
perpendicular to the glide cylinder) 


Ai 
© (V/200°6?) = Rae vaexDl = zm)" Of) [Sixes (10) 
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The left-hand side of (10) is the rate at which the number of vacancies 
forming the loop decreases (volume per vacancy =6*/+/2). The rate of 
emission of vacancies from the loop is then 
ay = 1/ 2rzvan eX - 7) moe he) 
di potas WO e/E ache a 
and this is independent of loop size. It is seen therefore that at a fixed. 


temperature the rate of change of area of a loop is constant. 


Fig. 10 


14 
2-4x10 


8x10- 


Histogram of loop sizes for the micrograph of fig. 1 of Hirsch, Sileox, Smallman 
and Westmacott (1958) taken from aluminium quenched from ~600°c. 
The average loop radius is 138A. The ordinate is loop density in 
number cm7%, assuming a foil thickness of 1000 A. See text for discussion 
of the radius 7) derived from Federighi’s results. 


Attempts to verify the form of eqn. (7) are inaccurate for the following 
reasons. The dislocation core itself is not imaged in the electron micro- 
scope but its strain field gives rise to a diffraction image lying to one side 
of the core (Hirsch, Howie and Whelan 1958, 1959). There is therefore 
an inherent inaccuracy in the measurement of loop diameters owing to the 
uncertainty in position of the dislocation. In addition, during heating, 
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the specimen may tilt, thus altering the contrast. The observations 
reported in figs. 3 and 4 are therefore considered to be in satisfactory 
agreement with theory. 


4.4, Annealing of Loops in the Bulk Material 


So far the chemical stress due tothe vacancy supersaturation produced 
by vacancies released from loops has been neglected. In, general it is not 
permissible to neglect this effect in experiments with the bulk material as 
the following considerations show. First consider the number of vacan- 
cies in equilibrium at various temperatures 7’ in a cube of aluminium 
1000 4 edge; this is approximately the volume occupied by a single loop 
in a region of typical loop density. Using the formula of Bradshaw and 
Pearson (1957) for vacancy concentration the following figures are obtained. 


Equilibrium number of 
vacancies per cube 


4 ope 1000 A edge 
600 2 LO= 
200 4 
170 id 

20 4x 107° 


The equilibrium number of vacancies near each loop in the annealing 
range (170—200°c) is therefore of-order unity, so that even if only a few 
vacancies per loop remain in the vicinity of the loop (or of other loops) 
appreciable supersaturation occurs. The degree of supersaturation, will 
be determined by several factors including the rate at which vacancies are 
produced by loops, the rate at which they move, and the nature and distri- 
bution of vacancy sinks (isolated dislocation lines, free surfaces and grain 
boundaries). The detailed solution of this problem is complicated. How- 
ever some simple generalizations can be made. The number of jumps, n, 
made by a vacancy in time ¢ is (Cottrell 1957) 


E 
N=Zvat exp ( - ch) 
and the mean square distance moved, assuming a random walk process, is 
oom E 
Aa? = b2zvat exp G ay dae Ne eg Or, Mea eee U9 
Then the time taken for a vacancy to move to a sink at a distance / from the 


loop at which it originates is 


i~ (Pera) exp (F) Ga a ett) 


Equation (11) also shows that if supersaturation effects are neglected the 
time interval ts between emission of successive vacancies by a loop is on 


average En+tk 
ts~ (v/2nea)"texp (=B5=4). cE re ws) 
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E 
Then tilts ~~ 0/ 2rl?ab exp (— Zp) ae eee es (0 HS, 


Now for thin foils the surface is a very effective vacancy sink; / is then 
about the average distance of a loop from the surface, i.e. ~ 1000 A. In 
bulk material however a vacancy may have to move about 1 to 10» before 
it encounters a sink such as a grain boundary or dislocation line (ignoring 
other loops which may act as sinks). The following are typical values of 
t/t; for aluminium at various temperatures obtained from eqn. (15) 
assuming a~ 70, Hy=0-76ey. 


Tc 1=10004 l=1p 1=10p 
170 8-7 x 10-2 8:7 870 
200 3x 10-1 30 3000 
220 6-6 x 10-4 66 6600 


It can be seen that for thin foils (J~ 1000 A), a vacancy from a loop will 
on, average reach the surface before the next vacancy is emitted by the same 
loop. The figures therefore suggest that the vacancy supersaturation, will 
be small and that the chemical stress can be neglected in this case. The 
theory of the previous paragraph is therefore expected to hold for thin foils. 

For loops in the interior of grains in bulk material however (/~ 101), 
t,/ts is large and increases with increasing annealing temperature. In this 
case supersaturation effects are expected to be important because vacancies 
are formed at loops at a rate greater than the rate at which they move to 
sinks. The effect of a vacancy supersaturation may be illustrated by using 
eqns. (1) and (2) to write eqn. (3) in the form 


OT esi SPL Ee Robe ke 
di = 5) ZV90 exp ( zn) E =| . . ° . (16) 


where Ac=c—cy. The rate of change of radius of a loop will therefore 
depend on the local supersaturation of vacancies and it is readily seen 
that a loop of radius less than re = a¢yb/ Ac will shrink while a loop of radius 
greater than re will grow. The large loops therefore grow at the expense 
of the small loops which disappear giving rise to a lower concentration of 
larger loops. This conclusion is in agreement with the observations on 
bulk material annealed at 190°c and 200°c (fig. 6) where a definite difference 
in the diameter and density of the loops in the as-quenched and annealed 
specimens was detected. 

In general local time-dependent inhomogeneities in vacancy concen- 
tration are to be expected. However the main physical features of 
the process as a whole can be illustrated in the following manner assuming 
a homogeneous supersaturation Ac. Consider first an idealized case of a 
metal containing no sinks for vacancies other than loops (e.g. a single 
crystal infinite in extent and perfect except for loops and vacancies). Let 
the metal contain a uniform density of loops of various sizes given by a 
distribution function g(r) such that there are g(r) dr loops per unit volume 
with radii between r and r+dr. Up-quench the metal to a temperature 
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: at 
T ~200°c. The rate of emission of vacancies now depends on loop size 
and vacancy concentration c. Since there are no vacancy sinks the 
following continuity equation holds: 


E 2 | Ac r | de 
Y ‘< oy ae ee al Fe 
4/ narnexp ( a). g(r)} « i ve; dr = sr (17) 


Now for small supersaturations once the initial ¢ has been established, 
both ¢ and de/dt are expected to be small in comparison with the two 
integrals subtracted on the left hand side of eqn. (17) at temperatures 
~170°c and 200°c, as substitution of typical figures suggests. With this 
assumption the right-hand side of eqn. (17) can be replaced by zero, giving 


ab | g(r) dr 
—~* ___ =ob/f . . . . . . (18) 
| rg (r) dr 

0 
where 7 is the mean radius of the distribution g(r). Electron micro- 
graphs suggest that 7 ~ 140 A in typical regions containing loops (see fig. 10). 
Putting «~70, ¢c/cy is found to be ~ 2-4 which confirms the assumption 
that c is small under near steady state conditions. For comparison 
the value of c/cy produced by quenching the metal from 600°c to room 
temperature is ~ 10°. 

Substitution of (18) in eqn. (16) gives 


Ona os 4 ee) (Ia 
di = Szvab axexp ( z)| 5 ; | e . . : (19) 


Equation (19) shows that the rate of shrinking of a loop now depends on 
r—f,1.e.%¢=F. Loops of above average radius will grow at the expense of 
those below average radius which will gradually disappear. Thus g(r) will 
change with time in such a way that 7 increases, i.e. Ac/cy decreases. 
Ultimately only a few very large loops will remain. 

In this idealized case the effects of grain boundaries and dislocations as 
vacancy sinks have been neglected. If these are taken into account even 
these large loops would shrink and vanish, provided they are not ulti- 
mately incorporated in a dislocation network. While therefore the final 
state produced by annealing bulk material for a sufficiently long time 
should be the same as that of the annealed thin foils (i.e. no loops), the 
detailed behaviour of loops as they anneal out will be different and is 
expected to depend markedly on their proximity to vacancy sinks. For 
example loops near grain boundaries would be expected to anneal out as 
in, thin foils, while those in the interior of grains should behave as outlined. 
In this last case since large loops are involved the annealing out time is 


expected to be longer. 


Acic,= 


§5. Discussion 
The agreement between observations and the theory of the previous 
section, leaves little doubt that the disappearance of loops is due to climb. 
The time and temperature dependences of the annealing rates in thin foils 


P.M. z, 
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are as predicted theoretically and are consistent with an activation energy 
of self diffusion of about 1:3ev. Bulk material experiments have shown 
results expected from the growth of large loops. 

It is of interest to consider the results in relation to the resistivity 
experiments of Panseri and Federighi (1958) and of Federighi (1959) on 
quenched aluminium. The first stage of their isochronal annealing 
curves (see fig. 2 of Panseri and Federighi) can be interpreted in terms of 
the removal of vacancies at dislocation loops (and possibly at small clusters 
of unspecified shape), and at fixed sinks such as network dislocations and 
grain boundaries. Now the temperature range in which the loops are 
observed to disappear is exactly the same as that found for stage Il 
annealing by Panseri and Federighi. Thus stage I] includes the annealing 
of resistivity due to the dislocation loops. It is therefore tempting to 
identify stage II with the annealing out of loops. It is not known with 
certainty at present whether very small clusters (other than loops) would 
be detected on transmission electron micrographs. No effects which can 
be attributed to very small clusters are observed. Even loops ~50A 
diameter would be expected to produce contrast in, the form of a dark spot 
about equal in diameter to the theoretical width of the image of a dis- 
location (~ 100A), and such dark spots have been observed (see fig. 1 of 
Hirsch, Sileox, Smallman and Westmacott 1958). However there is no 
evidence that these spots are due to anything other than loops. With 
these facts in mind stage II annealing in the quenched aluminium is 
tentatively considered to be due to the observed loops. Panseri and 
Federighi (1958) and Federighi (1959) have shown that the decrease in 
resistivity of aluminium during stage II annealing is ~3mpQcm after 
quenching from 600°c. If this is produced by loops of average radius 
~ 1404 and density ~ 101°cm~, the average resistivity of dislocations in 
loops is about 3 x 10-1? Oem per centimetre length of dislocation. It is 
emphasized that this value is an order of magnitude estimate only. To 
make this type of estimate worthwhile electron microscope observations 
and resistivity measurements should be made on the same specimens. 
However it is worth noting that the value is comparable with estimates of 
dislocation resistivities in other metals. For example Seeger (1957 ) 
quotes figures leading to dislocation resistivities of ~2:5x 10-9 and 
~ 10-18 A cm per centimetre length for copper and nickel respectively. 

Federighi (1959) has also studied the shape of the isothermal annealing 
curve of stage I resistivity in aluminium. He found (see fig. 4 of 
Federighi) that the fraction of resistivity remaining after annealing for 


time ¢ at temperature 7’ could be fitted very well for all except large t by 
the function. 


where mx 0:55 and 


E 
T,=T,9 OXp @ 
with B= 1-32ev. ois kT 
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T, represents a decay time ; more than 90% of the resistivity anneals out 
in time 7,. 7,9 was found to be 2-6x10-sec (see corrigendum to 
Federighi 1959). 

Now if the resistivity of a loop is taken to be proportional to its peri- 
meter (i.e. to its radius), the total resistivity of a distribution of loops g(r, #) 
will be proportional to 


| - rg(r, t) dr =F(t)ni(t), 


0 


where 7(¢) is the average radius of the distribution and n(t) the total number 
of loops at time ¢. On this interpretation therefore the explanation of 
Federighi’s empirical law should depend on the way in which the average 
radius and the total number of loops in the distribution changes with time. 
g(r, ¢) will of course be determined by considerations such as those outlined 
in §4.4 In the idealized case discussed there it is possible to formulate a 
differential -integral equation for g(7, t) ; however the equation seems soluble 
only by numerical methods. In the absence of a complete solution it is 
only possible to speculate. It is thought that the initial rapid decrease in 
Federighi’s annealing curve may be due to the initial rapid shrinking and 
disappearance of loops smaller than average radius and the growth of 
large ones. The subsequent slow decay would then be due to the slower 
shrinking of the large loops as vacancies are removed at network dislo- 
cations and grain boundaries at sinks. At present this problem is still 
unsolved and several other points are incompletely understood. It is not 
known for example what the contribution to the resistivity will be from 
heavily jogged dislocations. Furthermore examination of many specimens 
has shown that loops are not uniformly distributed over large regions. 
Loops tend to occur with high density in certain grains (usually free from 
isolated dislocations) and sometimes in patches in those grains. Other 
grains may show tangled networks of dislocation lines with lower loop 
density, e.g. fig. 5. These dislocations are probably produced by quenching 
stresses and are presumably heavily jogged by absorption of vacancies. 
During annealing these dislocations are observed to climb also; the jog 
density is thereby expected to decrease. This effect probably also contri- 
butes to resistivity although its magnitude is uncertain. Finally it must 
be pointed out that for small loops the contribution to electrical resistivity 
may not be proportional to 7; for example it may be proportional to the 
loop area, i.e. for?. The shape of the isothermal annealing curve is there- 
fore not determined by any simple considerations, but will depend on the 
complex annealing behaviour of loops in the bulk material. 

However a plausible explanation of Federighi’s 7, and 7,) can be obtained 
by noting that 7, depends on temperature in the same way as the theoretical 
quantity 7 previously defined (eqns. (7), (8) and (9)), which refers to a loop 
of radius r, in a thin foil. Equating 7, and Federighi’s value of 7,9. it is 
found that r,=4044. This is larger than the average loop size in material 
quenched from 600°c, as can be seen from the typical histogram of 


B2 
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loop sizes in fig. 10. However this result is in accordance with the 
qualitative ideas of annealing of loops in bulk material, i.e. that loops of 
larger than average radius grow at the expense of smaller ones. The 
magnitude of the decay time of the resistivity is therefore expected to be 
determined by the simple theory of § 4.3 applied to a loop of larger than 
average radius. It must also be noted that 7, is the time required for a 
large fraction of the resistivity to anneal out. It should therefore 
correspond to the lifetime of a large loop. 


§ 6. CONCLUSIONS 


Direct evidence for climb of dislocation loops has been obtained and the 
details confirm the theoretical predictions of the nature of the process, 
namely : 


(a) The loops disappear by emitting vacancies. 


(b) In thin foils the loops anneal out at an observable rate in the range 
of temperature > ~170°c and at a rate consistent with an activation 
energy of about 1-3ev. The rate has been followed and within experi- 
mental error is found to be as predicted by theory. Supersaturation 
effects are not important in thin foils owing to the very effective vacancy 
sinks provided by the surfaces. 

(c) In bulk material experiments the loops anneal out in the same 
temperature range, but the details are dependent on vacancy super- 
saturation effects. Large loops are found to grow at the expense of smaller 
ones. 

(d) The work has shown, that the loops disappear in the temperature 
range in, which certain electrical resistivity effects produced by quenching 
also anneal out. The loops are therefore tentatively identified as the main 
cause of the increase in resistivity, although at present other effects such as 
small clusters of vacancies (not observed in the electron microscope) and 
effects due to joggy dislocations cannot be completely ruled out. 


As regards the technique it appears that useful information can be 
obtained by high temperature examination of crystalline materials in the 
electron, microscope. Thus it has been possible to study directly for the 
first time the annealing kinetics of individual dislocations. Similar 
experiments on alloys may yield useful information on precipitation effects. 

Finally the fact that the annealing effects occur in about the same 
measured temperature range both in thin foils and in bulk material, is 
fairly convincing proof that under the conditions of operation used in this 
work the temperature rise of the thin foil due to heating by the electron 
beam is small (probably less than 15°c). This is also confirmed by the 
fact that loops usually disappear simultaneously both for regions inside 
and well away from the area exposed to the electron beam. 
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APPENDIX 
On THE QuESTION oF EMISSION VERSUS TRANSPORT OF VACANCIES AS 
THE RATE DETERMINING PRocESS OF THE CLIMB OF Loops IN THIN Foins 


Friedel (private communication) has pointed out to the authors that 
under certain, circumstances climb might be determined by transport of 
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vacancies by diffusion even in thin foils. He has considered two simple 
models—(a) diffusion of vacancies along a cylinder, (6) diffusion from the 
centre of a sphere. 
a) Diffusion Along a Cylinder 

Consider a loop of radius r at a distance / from both surfaces of the foil. 
Imagine the vacancies to be transported to the surfaces through cylinders 
of radius ry and length 7. Then the vacancy flux obtained from eqns. (1) 
and (2) (assuming c=¢, at the loop, i.e. no supersaturation) is 


Zvab7 exp ( — ak 
rs 2Qnrb dr ted © kT Real he F.b? 24 
by sack seaee TIES ic eck ED 


V2ar2 dt Vr 
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where V(=6?/4/2) is the vacancy volume, and where 
D=zvab? exp (— H/kT) 
is the self diffusion coefficient. However if emission occurs more rapidly 


than transport, c=c,exp (F'.b?/kT) at the loop and c=cy at the surface. 
The vacancy flux due to diffusion would then be 


D D Fb? 
= see | Pa ae 
Po Fes —— gradc= Ti [ exp (Sar TF ) ] (A 2) 


The climb rate is epee to be governed by transport of vacancies if 
> do, Le. if r/l<ce. Since ce~4, r would have to be comparable with 1 
for the climb to be determined by emission of vacancies. 


(b) Spherical Diffusion 

The above condition is however rather stringent and would be expected 
to apply to loops densely distributed in a plane at distance / from both 
surfaces. For randomly distributed loops in thin foils spherical diffusion, 
is probably a closer approximation. In this case a loop of radius r is 
imagined to be at the centre of a sphere of radius /(>r), i.e. for a very small 
loop in a thin foil. The vacancy flux at radius R for emission limited 
conditions (i.e. c=Cy at the loop) is now 


Dr Fb? 
$= spm’e| P(aGr)-} <a a SS ey 


However if c=cy exp (F.b2/k7) over a sphere of radius 7 at the loop and 
c=C, at R=1 (transport limited conditions) then 


e=ey-bes| exp ee )-. | 


and the vacancy flux due to diffusion is 


do= — ee | exp (ar) - 1. inde (A 4) 
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If the climb is determined by transport the condition ¢,>¢, gives Ce > 2. 
In this case therefore vacancy emission from the loops appears to be the 
rate determining process since ce~}. The actual case of loops in thin 
foils is expected to lie somewhere between these two cases but probably 
closer to case (6). Provided the loops are not too densely distributed and 
that the foil is not too thick, the approximation the authors have used in 
neglecting supersaturation effects in thin foils appears reasonable. 

Finally, it must be remarked that the total flux of vacancies from the 
neighbourhood of the loop (R~7r) under transport controlled conditions 
and assuming spherical symmetry is of the same order of magnitude as 
the total flux of vacancies emitted by the loop under emission, controlled 
conditions. The form of eqns. (5) and (6) therefore remains similar. 
Essentially similar conclusions have been reached by Weertman, (private 
communication) working on similar lines. It therefore appears that 
eqns. (5) and (6) are a valid description of the annealing process. However 
in the absence of more detailed observations on the rate of shrinking of 
loops it is not possible to decide by experiment whether the process is 
limited by emission or transport of vacancies. The considerations out- 
lined above and in §4.4 suggest that the rate would be determined by 
emission of vacancies. 
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ABSTRACT 


Dislocation loops formed by quenching thin aluminium specimens in 
iced brine were studied by the transmission electron microscope technique. 
As-quenched specimens contained loops 200 4 to 600 4 in diameter with a 
density of 101° loops per cubic centimetre. Quenched-in loops were 
completely removed by heating to 200°c for ten minutes or by cold rolling 
to a 5% reduction in thickness. The disappearance of the loops was 
accompanied by an increase in the number of irregular kinked dislocation 
lines. 


§ 1. [INTRODUCTION 


Tue equilibrium concentration of vacant lattice sites in a crystal decreases 
with decreasing temperature according to the approximate relation : 


ie 
C= exp (— =) 


where C,, is the fraction of the total number of lattice sites that are not 
occupied by an atom and U,, is the energy of formation of a vacancy. 
Therefore, when a metal is cooled, vacant lattice sites must be eliminated. 
There are only three ways excess vacancies can be removed: (1) migration 
to an external surface, (2) precipitation on an edge dislocation line, or (3) 
clustering to form small voids or prismatic dislocation loops as proposed 
by Seitz 1950). 

Resistivity measurements (Bauerle and Koehler 1957) suggest that 
when, the cooling rate is rapid enough, the concentration of vacancies present 
at some higher temperature can be frozenin. As long as the temperature 
is then kept low enough to prevent migration no change in vacancy 
distribution, can take place. 

When, an aluminium specimen which has been quenched from 600°o 
rapidly enough to retain the high temperature vacancy concentration is 
reheated to successively higher temperatures, the resistivity increment due 
to excess vacancies disappears in two distinct stages. The recent results 
of Panseri and Federighi (1958) are shown in fig. 1}. The first stage of 
recovery which accounts for 80° of the total resistivity increment lies in 
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+ Communicated by the Authors. 
{ All figures except fig. 1 are as plates. 
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the temperature range — 80°c to just below room temperature. As pointed 
out by Panseri and Federighi it is probably due to disappearance of 
individual vacancies by formation of clusters. Direct electron, microscope 
observations have shown that this thermal history does in fact produce 
vacancy clusters which collapse to form small dislocation loops 200 to 
600 A in diameter (Hirsch ef al. 1958). A dislocation density as high as 
10*°cm/cm$ can be produced in this way. Therefore the resistivity incre- 
ment remaining at room temperature may be due to the dislocations. — ‘The 
effect of the same thermal history on mechanical properties has been investi- 
gated by Maddin and Cottrell (1955). They observed a large increase in 
room temperature yield strength which must also be associated with the 
presence of prismatic dislocation loops. 


Fig. 1 
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Resistivity decay for isochronal annealing of 2 or 15 min at each temperature 
after quenching from 600°c. 


The second resistivity drop took place on heating into the temperature 
range between, 100°c and 200°c. The activation energy for the process 
was found to be 1:33ev (Federighi 1959). It seems likely that second- 
stage recovery may be associated with the disappearance of prismatic 
dislocation loops. Dislocation density could be greatly reduced by 
migration of vacant lattice sites from loops to other sinks or by growth of 
the largest loops at the expense of smaller ones. Small amounts of plastic 
deformation at room temperature also cause a decrease in resistivity 
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(Winterberger 1958). This suggests that loops can, be swept out by moving 
dislocations even at room temperature. The present experiments were 
undertaken to obtain, direct evidence concerning the mechanisms of this 
second stage resistivity decrease and of quench hardening. Changes in 
dislocation sub-structure in, quenched aluminium specimens caused by 
reheating to successively higher temperatures or resulting from various 
amounts of cold work after quenching were observed by electron 
transmission microscopy. 


§ 2. EXPERIMENTAL PROCEDURE 


Aluminium sheet specimens 0-005 cm thick of nominal purity 99-995% 
were quenched from 600°c into iced brine. The first series of specimens 
was then heated for ten minutes to either 100°c, 125°c, 150°C, 175°c, or 
200°c. The second series was cold rolled at room temperature to 24%, 
5%, 10%, 15%, or 45% reduction in thickness after quenching. A 
thicker specimen was cold rolled to 98% reduction in thickness. All 
specimens were then thinned by electropolishing in a solution of 20% 
perchloric acid 80% absolute alcohol to obtain fragments of about 1000 A 
to 15004 in thickness. The thin foil fragments were examined in an 
RCA-2E electron microscope operating at 50kv with either a 50 mil or a 
25 mil aperture. 


§ 3. REsuLTS AND DiIscuUSsSION 


As-quenched specimens which had not been heated above room 
temperature before electropolishing contained a dislocation sub-structure 
similar to that reported by Hirsch ef al. (1958). Figure 2 shows a typical 
field. Nearly all areas of the foils that had convenient thickness for 
observation, contained loops. They were absent only near occasional long 
dislocation, lines or within about a micron from a grain boundary. Loops 
appeared to be nearly circular and apparently did not enclose a stacking 
fault. Therefore the Burgers vectors were probably of the type (a/2)<110). 
The loops appeared to belong to four sets ; those in each family being quite 
accurately parallel to each other. This was particularly noticeable in 
foils oriented such that two sets were approximately at right angles to the 
plane of the foil because half of the loops then appeared as short straight 
lines. Foil orientation could not be determined because the instrument 
did not permit a diffraction pattern to be obtained. However, the observa- 
tions were consistent with the hypothesis that most loops lay approxi- 
mately on {111} planes. \ Neighbouring loops lying on a given (111) plane 
often were revealed with vastly different contrast. This might be expected 
since the loops on any one plane should have three different Burgers 
vectors e.g. (a/2)[110], (a/2)[101], (a/2)[011], on (111). A few appeared 
to have a dark centre and an approximately hexagonal shape. Whether 
this different appearance could be caused by a particular relationship 
between the Burgers vector, the crystallographic orientation of the foil, and 
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the direction of the electron beam, could not be determined because there 
was no provision for slightly tilting the specimen while under observation. 
An alternative explanation might be the presence of some loops enclosing 
a stacking fault. Most loops had the appearance shown at A in fig. 3. 
Hexagonal loops with uniform darkening across the entire area are shown 
at B. A few were revealed in double contrast asat C. The double contrast 
effect has been explained by Hirsch et al. (1960). 

Little or no change in the loops was produced by heating for ten minutes 
at 100°c. Figures 4, 5 and 6 show typical fields after heating to 125°c. 
As a result of a large number of observations, it was estimated that loops 
were absent in about one-third of the specimen volume. The number of 
irregular dislocation lines appeared to be greater and a few larger loops 
were found. Heating to 150°c decreased the proportion of specimen 
volume containing loops to about 50%. In the regions still containing 
loops their density was reduced and their average size increased compared 
to the as-quenched specimens as shown in fig. 7. Still fewer were present 
after heating to 175°c and essentially none remained after the 200°c heat 
treatment. Figure 8 shows a field still containing a few loops after 175°c 
heating. 

The results clearly showed that prismatic loops were eliminated by 
dislocation climb. Vacant lattice sites migrated from the smaller loops 
to the larger ones, causing growth of some at the expense of others or from 
loops to long dislocation lines. During this process there was a definite 
increase in the number of irregular kinked dislocation lines. Further 
heating caused disappearance of all loops and a straightening of the 
irregular kinked dislocations.’ The observed decrease in, dislocation, density 
took place in exactly the same temperature range as the second-stage 
resistivity recovery reported by Panseri and Federighi (1958). This 
supports the hypothesis that the decrease in dislocation density accounts 
for the resistivity decrease. The kinetics of the process is complicated by 
microscopic inhomogeneity of the specimens. Even in as-quenched 
specimens the density of loops apparently varied from one part of the 
specimen, to another. The disappearance of loops did not appear to take 
place primarily by a homogeneous coarsening process. Loops were 
eliminated in different parts of the same specimen at widely different 
temperatures. The higher the annealing temperature in the range between 
125°c and 200°c the greater was the percentage of specimen volume that 
was free of loops. In those parts of the specimen where loops remained 
longest, growth of the larger ones at the expense of the smaller ones did 
take place. The rapid disappearance of loops in some regions suggests 
that they may have been eliminated during the heating by a combination of 
glide and short distance climb of nearby long dislocation lines. By combi- 
nation of elements of the three dimensional dislocation network with 
vacancy loops the irregular kinked lines that were observed may have been 
produced. The regions that retained loops to the highest temperatures 
presumably were those farthest from network dislocations. 
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Further support for this general mechanism was provided by cold 
working after quenching. Dislocation loops combined at room tempera- 
tures with moving dislocations to form irregular kinked dislocation 
lines. All loops were swept out after only 5% reduction in thickness by 
cold rolling. A resistivity drop produced by room temperature defor- 
mation can therefore be associated with a decrease in total dislocation 
density. Figure 9 shows some loops remaining and irregular kinked 
dislocations after a plastic strain of 24% in tension. Figures 10 to 15 show 
the gradual change from tangles of kinked dislocation lines to a system of 
well-defined sub-grains as the amount of cold work was increased. 

The stress strain curve for quenched aluminium (Maddin and Cottrell 
1955) can, be qualitatively understood on the basis of these results. The 
yield stress was high, because the first dislocations to move appreciable 
distances became heavily jogged almost immediately by combination with 
loops. The resulting irregular dislocation lines were relatively immobile. 
At small strains the hardening effect due to formation of groups of tangled 
irregular dislocations was approximately compensated by the softening 
effect due to elimination, of vacancy loops; therefore there was little strain 
hardening. Dislocation density decreased during the first few percent 
plastic strain. At higher strains the rate of hardening and the changes in 
dislocation sub-structure appeared to be the same in quenched as in 
furnace-cooled specimens except that the hardness and dislocation 
structure reached at 5°% strain in the former was approximately equivalent 
to that produced by 10% strain in the latter. 


§ 4. CONCLUSIONS 


1, Dislocation loops formed by condensation of excess vacancies in 
quenched 99-995°% aluminium began to disappear on heating for 10min 
to 125°c and were completely eliminated by heating for 10min to 200°c. 

2. The process was not homogeneous throughout the specimen. Clear 
areas began to grow at about 125° and the proportion of specimen volume 
that was free of loops increased with increasing temperature of treatment. 

3. In the regions still containing loops after heating into the upper part 
of the range coarsening took place by growth of the larger loops at the 
expense of the smaller ones. 

4, Elimination of the high dislocation density associated with quenched- 
in loops probably explains the second-stage resistivity drop observed by 
Panseri and Federighi (1958), 

5. Irregular kinked dislocation lines appeared during the disappearance 
of loops and then, tended to straighten out and disappear again after the 
loops had been eliminated. 

6. Moving dislocations combined readily with the loops at room 
temperature. Irregular kinked dislocation lines were produced. A plastic 
deformation of 5% reduction in thickness by rolling was sufficient to sweep 
out all the loops—resulting in a decrease in total dislocation density. 
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ABSTRACT 


It is shown experimentally that nitrogen-locking is appreciably stronger 
than carbon-locking. This is contrary to the customary assumption and it is 
suggested that there may be chemical as well as elastic interaction with the 
dislocations. Boron is found to produce even stronger locking than nitrogen. 
These variations in the locking strength alter the ductile-brittle transi- 
tion temperature. The stronger the locking, the higher the transition 
temperature. 


§ 1. Locxrine By Boron 


Tue locking of dislocations in w-iron by carbon and nitrogen due to relief 
of elastic strain energy is now well known (Cottrell 1948, Cottrell and 
Bilby 1949). 

The present work began with an examination of the possibility of 
dislocation-locking by boron. Some doubt exists about the nature of 
the boron solution in «-iron. However, the observation of an internal 
friction peak at about room temperature and of a diffusion rate for boron 
in iron comparable to that of nitrogen and carbon suggests interstitial 
solution (Hasiguti and Kamoshita 1954, Thomas and Leak 1955, 
Samsonov and Tseitina 1955). Because of this and the fact that boron 
is a larger atom than either nitrogen or carbon, very strong dislocation- 
locking by boron seems possible. 

The approach adopted to obtain information about this possibility 
involved an examination, of the effect of boron on the lower yield point 
OLy.p, of a steel at 18°C and — 196°C over a wide range of grain diameter 1. 
This permitted a study of the effect of boron on the constants o, and k 
in the equation 

Sly p,=Ot+kl-?, a "es! te > wee (1) 


The form of this equation is consistent with the propagation of a Liiders 
band from grain, to grain at the lower yield point by the concentrated shear 
stress generated ahead of a slip band where it is held up by a grain boundary. 
The constant ’ is taken to be a measure of the shear stress tq required 


t+ Communicated by the Authors. 
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ahead of the slip band to unlock a dislocation. In the expression given 
in the Appendix, 


k~ 47'2aa, ees eat (2) 


) where r is the distance of the dislocation from the end of the slip band. 
The other constant op in (1) is the stress opposing the motion of an unlocked 
dislocation, along the slip plane. At room temperature, oy is largely due 
to the interaction, of the dislocations with fine precipitates, lattice defects 
and random impurity atoms, but at low temperatures another, very 
temperature-dependent contribution becomes dominant and this is 
apparently associated with the body-centred cubic structure (Cracknell 
and Petch 1955, Heslop and Petch 1956). 

To study the effect of boron, three steels were used, representing the 
addition of molybdenum and of boron and molybdenum to a 0-:13% 
carbon base material. The compositions are given in table 1 (A, B and 
C). This molybdenum—boron composition was chosen because it is 
representative of the principal boron steel used industrially. 


Table 1. Analysis, Weight % 


| re es _ “ Acid- © Acid- 
- we = 2 P Mo Al N 0 soluble’ B | insol.’ B 
A | 0-12 0-52 0-34 0-010 0-027 0-50 0:07 0-008 | 0-005 0-002 0-002 
rB | 0-13 | 0-53 | 0-40 | 0-011 0-030 0-48 0:07 | 0-008 | 0-004 <0-:0001 trace 
€ | 0-13 0-47 0-30 0-010 0-031 0-11 0-08 0-008 | 0-005 <0-0001 trace 
D | 0-12 0-51 0-02 0-050 0-029 a= <0-01 0-008 — <0-0001 nil 
-E | 0-019 | 0-005 | 0-006 | 0-014 | <0-002 | <0-01 | <0-005 | 0-007 —- nil nil 
-F | 0-045 | 0-47 0-02 0-035 0-018 — 0-01 0-006 — nil nil 


Ni, Cr, V, Cu, Co, Ti and W were all <0-02%, except for 0:08% Ni and 0:05% Cr in D. 


The general experimental procedure for the determination of the yield 
points followed that already described (Cracknell and Petch 1955). 
Specimens with 1 in. gauge length of 0-125 in. diameter were used and the 

_ strain rate was 2x 10-4sec-1. The specimens were machined, polished to 
3/0 emery and, to vary the grain size, they were annealed in vacuo in, the 
range 910°-1200°c, followed by cooling in the furnace at various rates. 
They were then tested. As commercially heat-treated, this boron steel 
is air-cooled and does not show a yield point probably because of the 
internal stresses arising from the low temperature of transformation from 
austenite. However, by using a cooling rate not exceeding 6°c min™, 
the yield point was retained. 

The results of these measurements are shown in figs. 1 and 2, and the 
values of o, and k obtained from them are in table 2. Throughout the 
present work, the least squares method was used. for curve-fitting. 

It is apparent that boron produces a substantial increase in the slope & 
both at 18° and —196°c. Thus, from (2), unless boron affects the 
distribution, of the dislocations so that r alters, this increase in k must 
mean that boron increases Tq, the dislocation-locking strength. 
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The lower yield point-grain size dependence at 18°c for steels A, B and C. 


G gives the limit of proportionality for the yield point-free material 
obtained from E. 


Table 2. o) and k Values 


Temperature G, k 
SMa °c 108 dynes em | 108 cgs units 
A (Mo-B) Seals 6-5 1-02 
B_ (Mo) 18 9-4 0-69 
C (killed) 18 7-4 0-65 
D_ (semi-killed) 18 el 0-74 
At (Mo-B) 18 6-6 0-66 
A (Mo-B) — 196 40-1 1-94 
B (Mo) — 196 45-0 1-28 
C (killed) —196 43-0 1-21 
D  (semi-killed) — 196 42-9 1-62 
EK’ (pure C) —196 49-3 0-98 
y (pure C) —196 44-5 1-01 
uw ( 


pure C-N) —196 44-5 1-59 


+ Annealed 750-850°c. 
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Table 2 also shows that when the grain size range in the molybdenum— 
boron steel A was established by annealing at 750-850°c, instead of above 
910°C, the effect of the boron on k& disappeared. The solubility of boron 
in «-iron falls rapidly below 910°c to <0-0004% at 710°o (McBride é¢ al. 
1954), so it seems that the amount of boron retained in solution in this 
annealing range was too small to affect the dislocation-locking. On 
cooling from the higher annealing temperatures, equilibrium was 
presumably not established and the boron was retained in solution. 
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The lower yield point-grain size dependence at —196°o for steels A, B and C. 
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§ 2. Locxine In KInLED AND SEMI-KILLED STEELS 
The steels used in §1 were fully killed (deoxidized) with Mn, Si and Al. 
Previously a semi-killed plain carbon steel had been examined (Cracknell 
and Petch 1955, Heslop and Petch 1956) and the o, and k values for this 
steel D are included in table 2. The composition is in table 1. It is clear 
that the slope k is appreciably greater for the semi-killed D than for,the 


fully-killed C, particularly at —196°c. 


Fig. 3 
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In an annealed semi-killed steel there is much more nitrogen than 
carbon in solution in the ferrite at room temperature and available for 
dislocation-locking. Fully-killing a steel, however, not only deoxidizes 
it, but in addition largely precipitates the nitrogen as AIN. Thus, this 
reduction in the slope & observed with the fully-killed steel suggests that 
removal of nitrogen from solution has weakened the dislocation-locking. 

To get some idea of the extent of the nitrogen removal in these steels, 
their strain-ageing was examined. This was carried out at 35°C, since the 
temperature is then too low for any ageing due to carbon. Specimens 
annealed to the same grain size were used. Figure 3 shows the measure- 
ments for the fully-killed C and the semi-killed D after 4% strain. Also 
included is a nitrogen-free, iron—carbon, alloy F’, described later. From 
the rate of ageing, it appears that the content of uncombined nitrogen in 
the killed steel C is about one-third of that in the semi-killed D. Figure 3 
also shows that the yield point elongation after ageing was shorter for the 
fully-killed steel. Since this elongation measures the strain required in 
one crystal to unlock a dislocation in the next, it seems that the nitrogen 
available in the fully-killed steel was insufficient for the development of 
full nitrogen-locking on ageing. 

A similar effect to that of full-killing upon & has been observed with 
high manganese additions (1—2°%) to a semi-killed steel (Heslop and Petch 
1957) and, although there is no precipitation of a nitride in that case, there 
is a 0-2ev attraction between the manganese and nitrogen atoms in 
solution, that may interfere with the supply of nitrogen atoms to the 
dislocations. 


§ 3. PurRE CARBON-LOCKING AND PuRE NITROGEN-LOCKING 


To examine further this effect of the availability of nitrogen on the 
dislocation-locking strength, pure carbon-locking and pure nitrogen-locking 
were prepared. ; 

Carbon and nitrogen were removed by wet hydrogen from specimens of 
a fairly pure Swedish iron and of a 0-47% Mn steel (E and F in table 1) 
until the yield point was suppressed}. Carbon alone was then replaced as 
described in§ 6. This gave two nitrogen-free steels, one with 0-20-0-26% C, 
the other with 0-09-0-12%C, 0:47% Mn. Their lower yield points at 
—196°c are shown in fig. 4 and the & values are in table 2 (E’ and EF’). 

It is clear that the pure carbon-locking gives a value of k only about 
_ two-thirds that of the semi-killed steel D at —196°c. This carbon & is 
even lower than & for the fully-killed steels B and C, which still contain 
‘some free nitrogen, and is close to the value reached in a high manganese 
(1:9%) steel (Heslop and Petch 1957). ae 

Replacement of the nitrogen to 0-008% in the carburized specimens 
returned & to the value for the semi-killed steel D (F” in fig. 4 and table 2). 
ee ee 

+ The limit of proportionality for the yield point-free iron obtained from E is 


included at G in fig. 1. 
02 
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Direct nitriding (0-013% N,) of the yield point-free material gave the same 
result with k=1-6 x 108 c.g.s. units at —196°c. (fig. 6). 

Thus, it seems well established that carbon-locking is weaker than 
nitrogen-locking and that the locking in an annealed semi -killed steel is 
predominantly due to its nitrogen. 


Fig. 4 
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The lower yield point-grain size dependence at —196°o after treatment to give 
pure carbon-locking (E’ and F”’) and after subsequent nitriding (F’). 


Specimens of F treated to give pure boron-locking (0-004—0-006% 
‘soluble’ B) have also been prepared (Codd, unpublished) and hore 
k value at 18°C is very close to that for the molybdenum-boron steel A 
so confirming that carbon-nitrogen-locking is replaced by boron-locking 
in this steel. , 
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§ 4. Toe REeLAtive STRENGTHS oF CarBon, NITROGEN AND 
BORON-LOCKING 


In (2), a change in the slope k can arise from a change in the locking 
strength rq or from a change in the distance r between the end of the slip 
band and the dislocation that is unlocked in the propagation of yield from 
one grain to the next. Replacement of nitrogen-locking by carbon-locking 
changes not only k but also its temperature-dependence. An alteration in 
would not have this effect. Also there seems no particular reason why r 
should change whenever the locking atoms change. On these grounds, 
the observed alterations in k have been interpreted as alterations in 
locking strength. 

In the original Cottrell—Bilby theory, the dislocation-locking in «-iron 
was thought to be due to a condensed atmosphere consisting of a single 
line of carbon or nitrogen atoms at the dislocation. It now appears that the 
atmosphere is more extensive than this (Cottrell 1953a, Dahl and Liicke 
1954, Thompson 1958). Thus, one possible explanation of the greater 
strength of nitrogen-locking compared to carbon-locking is simply that 
the higher solubility of nitrogen in a-iron below about 580°c may give rise 
to a higher concentration of impurity atoms at the dislocations. Against 
this idea, there is the consideration that escape from a stronger-locking, 
more concentrated atmosphere would be expected to require a greater 
activation energy and this would make nitrogen-locking less temperature- 
dependent than carbon-locking, whereas the opposite is observed. Also, 
when the carbon-locked steel F’ was quenched from 650°o, its k value 
remained unchanged at the low value (F”, fig. 5), and similarly on quenching 
the nitrogen-locked steel its k remained unchanged at the high value 
(fig. 6). Yet, in these treatments, the carbon concentration in the ferrite 
was greater than the nitrogen concentration. Thus, it does not seem 
possible to explain the difference between the carbon and nitrogen k values 
on the basis of a concentration effect. 

It is now suggested that the observed difference in the k values shows 
that the interaction energy between a carbon, atom and a dislocation is 
less than that between a nitrogen atom and a dislocation. Thisis, of course, 
contrary to the customary assumption based on the very nearly equal 
dilations produced by the two atoms and suggests that there must be some 
effect additional to the elastic interaction. 

Recent theoretical work has emphasized the possible directional character 
of metallic bonding (Altman et al. 1957) and the observations of what may 
be an appreciable Peierls-Nabarro stress in w-iron also suggests direction- 
ality of bonding (Heslop and Petch 1957). Possibly, therefore, in the 
locking of a dislocation in «-iron, the important effect additional to the 
EES SGT ae AO HITE. 9) SY RESO OE ALES By ELE SET ee es na ee 


+ There is an increase in the friction stress o) on quenching the iron—carbon 
alloy (fig. 5) because of the increased concentration of carbon retained in 
solution in the ferrite, but there is little effect on op for the iron—nitrogen alloy 
(fig. 6) because most of the nitrogen is already in solution in the annealed 
"specimens. 
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elastic interaction arises from a degree of residual chemical bonding at a 
dislocation to which the impurity atom becomes attached, as occurs in 
covalent silicon and germanium (Seitz 1952). 

There is some evidence that such chemical bonding might be stronger 
for nitrogen than for carbon. The free energy of formation of the iron 
carbides and nitrides are about the same (Quill 1950) in spite of the greater 
stability of the bonding in elementary nitrogen (Pauling 1948), which 
suggests that the iron—nitrogen, bond is stronger than the iron—carbon one. 


Fig. 5 
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The effect of quenching from 650°c (F”) on the strength of the pure 
carbon-locking at —196°c. 


Now consider the boron-locking. That it is the strongest is consistent 
with the atomic size of boron. The boron radius in «-iron is uncertain, but 
taking the value 0-944 from Fe,B, the tetragonal strain in a-iron should 
be ~0-53 and in the transverse directions the strain should be ~ 0-08 
Using the expression, due to Cochardt et al. (1955), these figures make the 
interaction energy with an atom situated directly below an edge dislocation 
~ 2-0 times as great for boron as for carbon. The observed ratio of & for 
boron and carbon is 2:0 at —196°c and ~ 1:7 at room temperature. This 
agreement is probably fortuitous since the boron radius is uncertain and 
chemical interaction is neglected. 


Locking of Dislocations in «-iron 39 


Fig. 6 
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The lower yield point—-grain size dependence at —196°c for pure 
nitrogen-locking; O annealed, @ quenched from 650°. 


§ 5. THe Friction STRESS AND THE DucTILE-BRITTLE TRANSITION 


The value of oy in (1) represents the friction opposing the movement of 
an unlocked dislocation along the slip plane. This is not a major concern 
in the present paper, so the observed values (table 2) will be discussed only 
briefly. 

The friction oy is due to the lattice itself and to interaction of the moving 
dislocation with precipitates, lattice defects and impurity atoms in solution. 
Between the semi-killed steel D and the fully-killed C there should be a fall 
in o, due to the removal of nitrogen from solution, but this appears to be 
compensated for by the increased interaction with silicon. The addition, 
of Mo in steel B makes a further contribution to o). The subsequent fall 
on adding boron (steel A) presumably arises from further nitrogen 
fixation and possibly from some interaction of the boron and molybdenum 
atoms. 

The variations in the dislocation -locking strength observed in the present 
measurements should be accompanied by variations in the ductile—brittle 
transition temperature (Heslop and Petch 1957, 1958, Cottrell 1958). 
The weaker the dislocation-locking, the easier is plastic deformation and 
the lower should be the transition temperature. This is confirmed for the 
present steels by fig. 7 in which the transition temperature is taken at 
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10 ft lb (~ 75% cleavage) for 5/16 in. V-notched specimens broken in @ 


Charpy-type impact test. 
The amount the transition temperature of the fully-killed C is below 


that of the semi-killed D agrees with previous measurements in which k 


was reduced by manganese (Heslop and Petch 1957). ; 
The 1-3 ton in-2 increase in o, between C and B (table 2) should give 
a 6°C rise in transition temperature (Heslop and Petch 1958). 


Fig. 7 
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The ductile-brittle transition temperatures for boron-locked (A), 
killed (B and C) and semi-killed (D) steels. 


§6. EXPRIMENTAL DETAILS 


For the preparation of the pure carbon-locking, tensile specimens 
(1 in. gauge length of 0-125 in. diameter) were partly machined from the 
Swedish iron E and from a similar material F with 0-5°% Mn added in 
the hope that it would help in the retention of a fine grain size. These 
specimens were heated in wet hydrogen at 850°c until samples showed 
that the carbon and nitrogen had been removed to such an extent that the 
yield point was suppressed (~ 600 hours). The treatment was continued 
for at least a further day and carbon only was then replaced by carburizing 
for a pre-determined time at 820°c in an atmosphere of purified hydrogen 
and hexane vapour. 

To remove nitrogen and other gases from the hexane, it was previously 
refluxed at reduced pressure (2mm) and temperature (0°c). After this 
treatment, the vacuum was released by purified hydrogen and, as an 
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additional precaution, the hydrogen was bubbled slowly through the hexane 
for some time before passing the hydrogen—hexane mixture into the 
carburizing tube. Since there was very little carbon pick-up until the 
hexane decomposed at ~815°c, the hydrogen-hexane mixture was used 
to protect the specimens during heating and cooling and the time at 820°c 
gave a reasonable estimate of the carbon ingress. After homogenizing 
for eight days at 1050°c in vacuo, the specimens were machined to final 
size, polished and then heat-treated in vacuo to establish the desired grain 
size. Each specimen was analysed for carbon. 

For the nitrogen-locking, nitriding was carried out at 580°C for 
~ 100 hours in a stream of hydrogen bubbled through 3% ammonia in 
distilled water. 


APPENDIX 


CALCULATION OF THE SLOPE k 


In Stroh’s (1955) calculation, the distance 7 between, the head of a slip 
band and the operated source (length s) in the next grain during the 
propagation of a Liiders band from one grain to the next at the lower yield 
point is taken to be at least equal to s. This is to allow sufficient room for 
the operation of the source. The shear stress ahead of a slip band depends 
upon 7—!/?, but the stress required for critical bowing of the source depends 
upon s~!, so, with s~r and no impurity-locking, the further the source is 
away from the slip band the longer it can be and the easier it is operated by 
the slip band. With a locked source, however, there will be a lower limit 74 
on the shear stress required for its working and the easiest source to bring 
into operation ahead of a slip band will then be the one for which the bowing 
stress and the locking stress are equal. 

Stroh thus obtains 

k ~ 2(27aub)"?, nes heheh, coe ogee ) 


where 1 is the rigidity modulus and b is the Burgers vector. 

This argument will however break down if the locking stress is greater 
than the bowing stress for all the sources present and it seems probable 
that this will be so in w-iron. At —196°0, tq is ~1 x 101°dynes cm~ and 
it is probably only a few times smaller at room temperature (Cottrell 
1958). Thus, the locking stress is ~ 100 times greater than that required 
for the critical bowing of sources of the length (10~4cm) normally expected 
in «-iron. ; 

The argument for (3) in a-iron has been supported by showing that it 
gives the correct temperature-dependence of k from the temperature- 
dependence of va. However, in this, ta was identified with the upper yield 
point, whereas it is now clear that the upper yield point is the stress 
required for the break-out of a plastic nucleus across the specimen 
(Cottrell 1953b) and not the stress for unlocking a source. The upper 
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yield point is influenced by the stress concentration ahead of slip bands 
in the plastic nucleus, so it depends on grain size in a similar way to the 
lower yield point (eqn. (1)). Thus the temperature-dependence of the 
upper yield point arises principally from the temperature-dependence of 
the friction stress o, and it does not directly give the temperature- 
dependence of the locking strength vq in the way assumed in this support 
of (3). 

oe of a value of r determined according to Stroh’s suggestion, it 
may be that r is simply the average distance between the end of a slip band 
and the nearest source in the next grain. Then, from the equivalence of 
a crack and a slip band in the relaxation of shear stress 


3(O1y.p, — %) (1/47)? = Ta, 
or 
k=4r'l27q. oi sikh (4) 
At —196°c, with ta~1x 10dynesem~ and k=1-6x 10% c.g.s. units, 
ris ~ 0-2 x 10-* cm. 

The main difficulty with (4) is that r has to be constant, independent 
of grain size. It may be that the width of the grain boundary stress field 
imposes such a constancy. 

Another possibility is that instead of operating a source in the next 
grain, the yield is propagated by forcing a length of dislocation out from 
the grain boundary. 
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ABSTRACT 


The mechanism of dislocation loop formation by vacancy condensation is 
examined in a quantitative manner. The critical vacancy supersaturation 
necessary for the growth of the loops and the climb rate of the dislocations 
are determined. It is shown that the mechanism proposed for lineage 
formation based upon loop growth (Frank, 1956a) is inapplicable. The 
number of dislocation loops produced by collapsing vacancy discs is calculated 
as a function of the cooling rate, cooling range and dislocation density 
existing in the crystal before cooling. 


§ 1. IntTRoDUCTION 


ESSENTIALLY five different mechanisms exist for the introduction of 
dislocations into freshly grown crystals: growth accidents (Frank 1950), 
vacancy condensation (Frank 1950, Seitz 1950, Kuhlmann-Wilsdorf 1958), 
segregation of solutes (Frank 1950, Friedel 1956, Goss e¢ al. 1956, Tiller 
1958), phase transformation and precipitation (Brooks 1952) and thermal 
and external stresses (Billig 1956, Wagner 1958). If one considers single 
phase crystals grown carefully from the melt so that both thermal and 
external stresses are absent we may neglect the last two mechanisms. In 
general, then, one needs to consider only the first three mechanisms in 
order to determine the total dislocation, concentration either grown into 
the crystal at the solid—liquid interface or quenched-in during cooling of 
the crystal. 

One would expect that not only the dislocation density but the 
dislocation arrangement in the crystal would differ depending on the mecha- 
nism by which they are introduced into the crystal. However, except — 
in a few cases, little quantitative theoretical work has been carried out on 
the genesis of the original dislocations in metal crystals. Most theoretical 
investigations have been concerned with the somewhat simpler problem of 
the stable arrangements of dislocations either in two dimensional arrays 
(boundaries) or in three-dimensional networks. ' 

Recently, Kuhlman-Wilsdorf (1958) has reconsidered the formation of 
dislocation rings from collapsed vacancy discs and proposed that these 


panemeeenan ering: se en Se Se 


Communicated by the Authors. 
| Now at CNEA, he. Lib. Gen San Martin 8250, Buenos Aires, Argentina. 


44 G. Schoeck and W. A. Tiller on 


dislocation rings are actually the main sources for the multiplication of 
dislocations during plastic deformation. From the kinetics of the 
electrical resistivity recovery and from quench hardening, it has been 
concluded by Kimura et al. (1959) that the dislocation rings can be formed 
from vacancy aggregation provided the crystal is quenched over a suffi- 
ciently large temperature range. More recently Hirsch et al. (1958) have 
actually observed such dislocation loops in quenched Al (Hirsch et al. 1958), 
and defects related to dislocation loops in the noble metals (Hirsch and 
Silcox 1958). 


Fig. 1 


Striation-type dislocation array (after Frank). 


The vacancy condensation model has also been used by Teghtsoonian 
and Chalmers (1951, 1952) and Atwater and Chalmers (1957) to account 
for a particular type of lineage structure called ‘striations’ observed to 
occur during the growth of crystals from the melt. This sub-structure 
consists of an assembly of low-angle grain boundaries which partitions the 
crystal into rod-like elements that are generally parallel to the crystal axis. 
The striations are often of the order of 1 sq. mm in cross-sectional area, the 
different rods being rotated with respect to each other about the growth 
axis by angles varying between a fraction of a degree and several degrees, 
as illustrated in fig. 1. In growing a single crystal without the use of 

a seed it is observed that the striations do not form immediately. They 
are observable only after a considerable length (approximately 1 cm) of 
solid has formed. After this incubation distance, the misorientations 
appear spontaneously and increase to a stable value with a little further 
solidification, 

Frank (1956 a) has developed a detailed picture of the vacancy collapse 
mechanism to account for this ‘striation’ sub-structure giving a quali- 
tative explanation of the origin, size and initial misorientation. It is 
very easy to show that the dislocations in the lineage array could not have 
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been formed by the collapse of vacancy discs since the total vacancy 
concentration of the crystal is 2 to 3 orders of magnitude too small to 
provide the lineage size and misorientation. Therefore, Frank postulated 
the following. 

The edge dislocation loops, formed by vacancy condensation at some 
distance behind the interface, grow in size to reduce the vacancy super- 
saturation. Those loops which lie in planes parallel to the growth axis 
will grow towards the interface and catch up to it. Once the dislocation 
has entered the interface it will contribute to the misorientation of the 
growing solid without further vacancy condensation. This is illustrated 
in fig. 2. He maintained that the elastic interaction between edge 
dislocations should cause them to assemble into parallel arrays as shown. 
A system of these arrays leads to the lineage type substructure (fig. 1). 
The ends of the loops in fig. 2 (a) are in an unstable position and are attracted 
into the interface with the liquid (fig. 2(b)). The dislocations are then 


Fig. 2 
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(a) (b) (c) 


The climbing of dislocation loops into the solid-liquid interface. 


no longer effective for removing vacancies merely by lengthening and the 
vacancy excess will drive them apart. A new dislocation loop will then 
form preferentially in this central region (fig. 2 (c)) and the amount of 
rotation of the lineage boundaries will continually increase until the rate 
of production of new loops is balanced by the rate of annihilation of 
dislocations by the meeting of opposite pairs. 

The purpose of the present paper is to examine the vacancy conden - 
sation mechanism in a quantitative manner tv determine (i) the critical 


46 G. Schoeck and W. A. Tiller on 


vacancy supersaturation necessary for the formation and growth of 
dislocation, loops by vacancy condensation and (ii) the climb rate of the 
dislocations as a function of the vacancy supersaturation. From these 
calculations and the experimental observations on, striations it will be 
shown, that Frank’s model cannot account for the sub-structure. From 
the calculations it is also possible to determine the quenching rate and 
quenching range necessary for the formation of the dislocation loops by 
vacancy condensation. Finally, it is possible to estimate the density of 
dislocations introduced into the crystal by this mechanism during crystal 
growth. 


§ 2. CONDENSATION OF VACANCIES 


If the concentration of vacancies in a solid exceeds the thermal 
equilibrium value there is a tendency for the vacancies to anneal out by 
condensing into larger aggregates. The determination of the vacancy 
super-saturation necessary to nucleate this process is a very complex 
problem and has not yet been treated. It is very likely that the nucleation 
takes place heterogeneously in the neighbourhood of solute or impurity 
atoms rather than homogeneously, but it is not possible to state whether 
these aggregates will be preferentially pores or flat discs (Frank 1956 b). 
However, since we are interested in the formation of dislocations, we may 
assume that vacancy discs are formed which collapse into dislocation rings 
after they reach a critical size. After this has happened we can calculate 
the supersaturation necessary to keep the dislocation rings growing to a 
larger size by climbing. We thus know that the actual supersaturation 
must be at least this value. Since all dislocations formed by vacancy 
condensation must go through this early stage of growth we will obtain, in 
this way, a lower limit for the necessary vacancy supersaturation. 

As has been, discussed by various authors (Kuhlmann-Wilsdorf 1958, 
Kimura et al. 1959) the collapse of the disc will take place when the energy of 
the disc (which is essentially the surface energy) exceeds the energy of the 
dislocation ring (which is essentially strain energy). Usually the collapse 
of the disc will also form a stacking fault since the two lattice planes 
brought into contact with each other will generally have an incorrect 
stacking sequence. In very few cases the surface energy of this stacking 
fault is known. Since a coherency exists across the interface, the stacking- 
fault energy should be of the order of a coherent twin-boundary energy 
and it will not exceed, and in many cases be considerably less than, the 
energy of a large angle grain boundary. It seems therefore safe to 
assume that the stacking-fault energy yst is smaller than about 4 or 4 of the 
surface energy y. To obtain an upper limit to the collapse radius in our 
calculations, we assume for numerical calculations ys¢=4y in all cases where 
Yst iS unknown. 

If r is the radius of the vacancy disc, the energy of the dise Hq for not 
too small values of 7, will be 


Eax 2nry, ° e ° ° e ° =). ° e (1) 
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The energy of a corresponding dislocation ring with the Burgers vector b 
perpendicular to the plane of the ring (Franz and Kroener 1955) and 
enclosing a stacking fault of surface energy ys will be 


1 Sr 2 
B= F — by Ee ~ i + Qa | 733 | + Tryst, sore G,) 
where G is the rigidity modulus, v is Poisson’s ratio, fy is the radius of the 
dislocation core-and L is the average strain energy per atom of this core. 
For numerical purposes we use v=}. 

We may determine the minimum core radius by noting that the free 
energy of the core cannot be greater than the free energy of a cylinder of 
liquid of radius Ry. If one determines the core radius using such an 
energy balance including a term which incorporates the solid—liquid inter- 
facial energy it is found that Ry~b/2. However, for such a small cylinder 
of liquid the surface energy term should probably be excluded. If this 
term is neglected the core radius is given by 


is el GIAH U2 (Tm) ¥ es 
b = 2x(1—v) SAT | SAI Gyo tty or mils? 


where AZ is the enthalpy difference between liquid and solid, 7m is the 
melting temperature and AJ =7,—7. For a typical value of 


AT /T m~ 0:25, K/b=2 
and eqn. (2a) can be put in the form 


tees r 
i= Bresette | me +2] ALaTi ata e eth Seat, P20) 
For &,/b=2, 1 and 0-5, Z=1-98, 1-48 and 1-96 respectively. We will 
choose the value of Z=1-8 in our calculations as a representative value. 
The condition for collapsing, H,= Eq, gives the critical radius r, from 


the equation 
b r o |. 42(2v— vst) 3 
a | toe) +1 3 = sreygeer wet eee Os (3) 


The critical radius 7, for several metals is given in table 1. In several 
cases 7 is astonishingly small and collapsing would occur according to 
eqn. (3) if the diameter of the disc is of the order of a few atomic distances. 
This estimate may not be quite realistic since the equation for the 
dislocation energy and the concept of surface energy may lose their signifi- 
cance for such small dimensions. Nevertheless it seems safe to assume 
that a disc with a radius of 10 atomic distances will have collapsed to form 
a dislocation ring. For some numerical estimates we will later assume 
that a vacancy disc exceeding a radius of 10 or 15 atomic distances will 
have formed a dislocation ring. Once the dislocation ring is formed, it is 
under the action of its own line tension and can grow or shrink by climbing. 
It will expand only if there is a sufficient supersaturation of vacancies. A 
vacancy condensing at the dislocation ring will increase its radius and 
hence its strain energy H, At the same time the free energy Fy of the 
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crystal due to the vacancy supersaturation, will decrease. If n is the 
number of vacancies, the critical supersaturation necessary to keep a 
dislocation ring of radius 7 from shrinking is obtained from the condition 


OL. on 
—— <0 tb eee err ee eee ee) 
ar Yar ~ ( 
where ty is the chemical potential of a vacancy which is given by 
ory Cc 
eae ee ET Se ol eee Te a 
sae ore ne, (6) 


Here c¢ is the actual vacancy concentration and c) the concentration of 
vacancies in thermal equilibrium at temperature 7’. 


Table 1 


G b y 
1011 dyne/em? | 10-8 cm | erg/cm? | erg/cm? 


6257 <ty 


t+ U.S. Atomic Energy Comm. Publ., (Liquid Metals Handbook, Wash. 1952). 
t J.C. Fisher and C. G. Dunn, Imperfections in Nearly Perfect Crystals (John 
Wiley & Sons, Inc., New York 1950). 


Since we want to obtain only a lower limit for the necessary super-.- 
saturation we are justified in neglecting the term in #, in eqn. (2) which 
contains the energy yst of the stacking fault when we evaluate eqn. (5). 
The result then becomes somewhat simpler and we obtain from eqns. (2c), 


(5) and (6) ‘ 
ne = oye | B(E) +28 J. ere) 


In order to obtain an estimate of the temperature at which the 
dislocation ring will be stable we must know the actual concentration c 
of vacancies present. The concentration of vacancies in thermal 
equilibrium at temperature 7’ is 


Cp= Boxp (=. Ua/kT)s. 15a eee (O) 
where U; is the energy of formation of a vacancy and the constant B which 
is of the order of unity may contain an entropy term. Let us now assume 
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that the equilibrium concentration of vacancies at the melting point, 7m, 
is frozen in at the advancing liquid—solid interface. Since the mobility of 
vacancies at 7’ is extremely large it is very unlikely that any in excess of 
this number are frozen-in. Only if no sinks are present will this number 
be retained in the interior of the solid during cooling to a temperature 7 
and hence we obtain as an upper limit for c 


e—Bexp{— Us/kT an) i se ee eS) 

If we call AJ’ =7,—7' we obtain finally from eqn. (7) with eqns. (8) and 
(9) 

AVE SI acre +f 
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(The entropy term 8 drops out if it is temperature independent.) Equation 
(10) gives the relative undercooling A7'/7Tm necessary to keep a dislocation 
ring of radius r from shrinking. The magnitude of A7'/7' depends on the 
parameter a=3Gb?/87U;, which is listed in table 2 for various metals. 
When measurements of U; were available the correct value was used 
whereas otherwise the value of U; was assumed to be half the activation 
energy for self-diffusion. Equation (10) is plotted in fig. 3as AT'/7'p vs. r/b 
for various values of x. The resulting supercooling for a dislocation loop 
radius of r= 10b and r= 15d is then also listed in table 2. 

To obtain an estimate of the order of magnitude for A7'/7'm we may note 
that U; is of the order (Friedel 1956 b) 


UO. bolero t meet (1a 


Using this estimate we obtain as a typical value, 1a=0-6. For this value of 
a in fig. 3 we see that in order to keep a dislocation of radius r= 106 from 
shrinking, the relative temperature drop would have to be A7'/T'm 0-35. 
The corresponding value for a dislocation ring of radius r= 15b would be 
AT /Tm~ 0-21. Assuming that a radius of r=100 is realistic for a newly 
formed dislocation ring, the relative supercooling varies between 16% and 
50% for the various metals (table 2). This means that the temperature 
at which these dislocation rings are stable is about 100°c to 400°C below 
the melting temperature if wé are considering crystal growth or the 
temperature from which they are quenched if we are considering a quench- 
ing experiment. It must be remembered that these calculations assume 
the existence of no other vacancy sinks in the solid, otherwise the vacancy 


+ It has been pointed out by the referee that if the vacancies condense to 
form a ring of Frank sessile dislocation including a stacking fault, the Burgers 
vector of whole dislocations, b, should not be used. Equation (10) and the 
tabulated values of A7'/T' » refer, therefore, to the case of no stacking fault. 
However, if one wishes to treat the stacking-fault case, the right-hand side of 
eqn. (10) must be multiplied by two-thirds to take account of the reduced 
Burgers vector and a term b*yst/Us added to take account of the stacking-fault 
energy. For rings of r=10b—15b the two changes tend to cancel for reasonable 
yst leaving AT'/T'm negligibly changed. 
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supersaturation at any temperature will be lowered and the (A7/7'm) 
required to keep a loop from shrinking must be increased. 


Fig. 3 


l2 14 16 ISit 520 


lO 

r/b 

The relative undercooling, A7'/7’m, necessary to keep a dislocation ring of 
radius r atoms from shrinking (7, =melting temperature). 


§ 3. RATE oF DIsLocATION CLIMBING 


It was noted earlier that Frank’s mechanism for striation formation 
(Frank 1956 a) requires that the dislocation loops, once formed behind the 
interface, must be able to climb towards the moving interface at such a 
rate that they can catch up with it. Therefore, we want now to consider 
the climb rate of the dislocations. 

To simplify the problem we shall consider a crystal at a constant 
temperature, 7’, with a constant excess concentration, ¢—¢», of vacancies 
present. In the crystal there may be a straight dislocation line which 
forms a perfect sink for vacancies. This last condition may perhaps not be 
satisfied for dislocations with high jog energies and hence a low concen- 
tration, of thermal jogs. However, if the diffusion rate along the dislo- 
cation in, its immediate neighbourhood is greatly enhanced, even a small 
number of jogs will be able to maintain the concentration of vacancies in 
thermal equilibrium all along the dislocation line. Furthermore, we are 
mainly interested in obtaining the fastest rate a dislocation could possibly 
climb ; thus, we may neglect the problem of jog formation which could only 
slow down the climb process. We assume therefore that the concentration, 
of vacancies at the dislocation will be kept in thermal equilibrium, cy, and 
that the excess concentration will anneal out there. There are two reasons 
why vacancies diffuse to the dislocation, line: (i) there is a chemical force 
due to the concentration gradient directed towards the line sink and (i1) 
the dislocation line attracts vacancies due to an elastic interaction (Friedel 
1956b). However, since the elastic interaction. energy is proportional to 
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+, the attractive force falls off rapidly with distance. As will be shown 
in, Appendix I, in most cases the elastic interaction is too small to effect 
the steady-state climb rate. 

From the formulation, of the physical problem it is obvious that the 
dislocation will reach a steady state where it climbs with constant velocity. 
If we neglect the elastic attractive force an exact solution, can be obtained 
(Jacob 1956) which is given in Appendix I]. In order to gain a better 
understanding of how seriously our simplifying assumptions affect the 
final result we will derive the magnitude of this velocity by a simple physical 
argument, the result differing from the exact solution only by a numerical 
factor. 

The dislocation will drain the surrounding region from its excess- 
vacancy concentration, c—cy. After a radius ry has been drained, the 
number m of vacancies having arrived at the dislocation per atomic 
distance 6 are 


=n (C= 6). | see 


After m vacancies have arrived per atomic distance 6 on the dislocation 
line the distance it has climbed will be 


t=mbi-+) .« 29s 6 eee) 


If /=r, the dislocation will essentially have climbed out of the region 
exhausted of vacancies and moved into a region, which has not yet been 
drained. Hence from eqns. (12) and (13) we obtain as the radius of the 
drained region 

b 


m1(C—Cp) ; i) 


foe 

The steady state velocity of climbing vy will now essentially be determined 
by the time fy it takes to drain an area of radius 7 and will be of the order 
Of Vg=7o/to. Lf the dislocation tends to move much faster insufficient 
vacancies diffuse to it to maintain the higher rate, and if it tends to move 
much slower the vacancies which arrive would push it ahead faster. Now, 
in the time ¢ a vacancy can move over a distance x=(Dyt)"?, where Dy is 
the diffusion coefficient for vacancies. Hence the radius 7) which is 
effectively drained of vacancies will be smaller but of the order of magnitude 
19 * (Dyty)'”. By use of eqn. (14) we obtained then 
%% _ TDy(c— Cp) 
roe b ea 
‘The exact solution in Appendix IT shows that the right side of eqn. (15) 
must be multiplied by a numerical factor of f=0-2—0-1. If we further 


notice that the coefficient for self-diffusion, is given by D= CoDy (assuming 
@ vacancy mechanism for self-diffusion) we obtain finally 
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We must keep in mind that eqn. (16) has been derived for straight 
dislocation lines. It will therefore give a reliable result only if the radius 
of curvature of the dislocation line is larger than the radius ry of the drained 
region of eqn. (14). It can be shown that the climbing rate of a dislocation 
ring with radius 7 much smaller than 7, will be smaller than Vo if we neglect 
the action of the attractive potential. However, for much smaller r, the 
attraction of vacancies to the dislocation line becomes important, but we 
would not expect v, to be very different from that already calculated. An 
exact solution for an expanding ring would be very complicated and no 
attempt has been made to calculate the climbing rate in this case. 

If we use eqn. (16), we see that a dislocation ring of 10 atoms radius is 
generally climbing, at its equilibrium temperature, towards the interface at 
a rate slower than a typical rate of advance of the interface (10-2-10-3 cm / 
sec). However, a dislocation of 15 atoms radius is climbing, at its equili- 
brium temperature, at a rate which is equivalent to the rate of advance 
of the interface. 


§ 4. FORMATION OF STRIATIONS 


1. A critical test of the origin of striations by vacancy collapse is to 
grow crystals under an ambient temperature so high that the critical super- 
saturation for the growth of dislocation rings is never reached. In this 
laboratory, crystals of Pb and Sn were grown in a horizontal crucible in an 
ambient temperature maintained at 20°c below the melting point. These 
crystals exhibited a fully developed striation type sub-structure. This. 
sub-structure could be seen, on the outer surfaces of etched crystals and on 
etched sections perpendicular to the growth axis. Since in these experi- 
ments A7T'/T 20-04, the radius of a ring which can expand by climbing 
would have to be at least 35b in Pb and 85bin Sn. It seems safe to say 
that rings of this size cannot be formed by spontaneous nucleation. Hence, 
we may conclude that collapsing vacancy discs are not responsible for 
striation formation. Mie 

Other objections to Frank’s model (1956 a) can be made by considering 
the experiments of Teghtsoonian and Chalmers (1951, 1952). From the 
calculations of AT/7'm necessary for a dislocation ring of 10 or 15 atoms. 
radius to grow, we can, see that the loops must nucleate several centimetres 
behind the solid—liquid interface. If we denote this value of A7/7Tm by 
(AT/Tm)c and the temperature gradient in the solid as Gs(°C/ em), then the 
distance, L, behind the interface at which the loops nucleate is given by 
Ls (AT/Tm)c(T'm/Gs). If we assume, for the moment, that Frank s 
(1956 a) model is correct and that all the vacancies collapse to form dis- 
locations in lineage arrays, the angular misorientation, 0, between the 
striations at Lem from the beginning of the crystal (L is the SCs 
distance), is given by tg9=c. Thus, for reasonable at of Gs=15°c ae 
and 30°c/em for Pb and Cu respectively and for cx 10™ we see that the 
incubation distance for striation formation should be about 5-10 cm with 
a misorientation of about 10-2degrees. However, the incubation distance 
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is observed to be less than about lem at high growth rates and the 
misorientation is observed to be about 1 degree. 

2. We should now consider if a loop, once formed in the crystal, can 
catch up with the moving interface. As we see from table 2, a dislocation 
ring of 10 atoms radius will climb, at its equilibrium temperature, at a rate 
slower than the rate of advance of the interface and thus will be left behind. 

However, a dislocation of 15 atoms radius is climbing at a rate which 
is equivalent to the rate of advance of the interface. If the initial loop 
size is 15 atoms radius there is some possibility that the loops will grow to 
a higher temperature if the interface is moving slowly. Whether it will 
climb close enough to the interface to be pulled into it by image forces 
must be decided. 

If we write the diffusion coefficient in the form D =D, exp (— Ua/Rt) the 
climbing velocity v» of eqn. (16) can be put in the form 


‘5D U U 1 1 = 


and using the empirical relations Uq~397'm, U;x4Uaq (Birchendal 1951), 
the reduced velocity V =v,b/D, is plotted vs. Tm/T in fig. 4. We see that 
there is a maximum in V at AT'/T~0-1. The climbing velocity vanishes 
for small 7’ because there the vacancies become immobile and it also 
vanishes at the melting point, 7m, because there the supersaturation 
disappears. For A7'/Tmaz10- the velocity is practically linear with 
temperature. 

Before a dislocation line can experience a significant influence of the 
image force at the solid—liquid interface it must move to, at least, a distance 
of 10-?cm from the interface. Therefore, we are interested in v, at this 
point. With Gs=15°c/em and 30°c/em for Pb and Cu respectively, 
AT'/T m at this point is given by AT/Tm= 10? Gs/Tm ¥ 2:5 x 10-4 for both 
these metals. Thus, at this point behind the interface the climb rate is of 
the order of 10-*-10-> cmsec and hence is about two orders of magnitude 
smaller than the conventional freezing rates. One must then conclude 
that the dislocations cannot catch up to the interface but can only trail 
behind it at a distance such that vp equals the freezing rate. From fig. 4, 
we see that this occurs at A7'/7mx.10-* which corresponds to a distance 
of about 0-5cm from the interface. These calculations neglect any effect 
that thermal stresses may have upon the climb rate. 

3. We have shown by various arguments that striations cannot be formed 
from collapsed vacancy discs. However, it is still possible that a number 
of independent dislocation loops may be produced by climb towards the 
interface in a manner similar to that suggested by Frank (1956a). Let us 
consider this modified picture and calculate the upper limits to the number 
of dislocations introduced into the crystal in this manner. We will 
continue to assume that no other vacancy sinks are present in the crystal. 
Considering loops moving with the interface at a value of AT [Pay= 1073; 
from eqn. (10) the minimum radius of the leading edge of the loop, 7, may 
be determined to be about 10¢atoms. From eqns. (8) and (9), the vacancy 
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excess is 10-°; thus, the maximum number of elongated loops of width 


2r in the crystal at the position correspondi 
ading to AT/Tm=10-2 : 
given by : i (aa 
C—Cy 
Pmax = orb (18) 


Climb Rate > Freezing Rate 


a Climb Rate <Freezing Rate 


| | 1.2 eS \.4 PS 1.6 
Teal: 


The reduced velocity of climb, V, as a function of T'y,/T'(V =05b/D, where 
v)=climb rate, b=Burgers vector and D,=diffusion constant). 


which gives pmax © 10®cm~*. From this calculation we can see that, if no 
other vacancy sinks are present in the crystal and if the dislocation loops 
form at a value of (AT/T'm)c where the climb rate is greater than the rate 
of interface advance, then a relatively large number of dislocation loops 
- will form and trail the interface at a distance such that the climb rate is 
equal to the rate of interface advance. From Dash’s (1958) work on 
dislocation-free Si crystals it would appear that ( AT'/T'm)c is too large to 
allow the necessary climb rates. 


56 G. Schoeck and W. A. Tiller on 


4. It has been shown experimentally and theoretically that the “vacancy 
dise mechanism’ cannot account for the striation-type sub-structure, so 
the cause of striation formation is still unknown. It has been, suggested 
that they most likely form from the dislocation, lines grown into the crystal 
at microsegregation boundaries (Friedel 1956 a, Tiller 1958). The density 
of these dislocations is of the order of 10®-108 lines/em?, which is sufficient 
to produce the striation arrays. 

Evidence of these dislocation lines have been observed by etch pit 
techniques, in Zn (Bédcek and Kratochvil 1959, Damiano private com- 
munication). Aggregation of these dislocations into micro-arrays occurs 
behind the solid—liquid interface, the arrays close to the interface being 
an order of magnitude smaller than striations and increase in size as the 
distance from the interface increases (Damiano, private communication). 

A driving force for the growth of micro-arrays into striations is given by 
the strain energy of the system since the energy of the crystal containing a 
fixed number of dislocations arranged in small-angle cell boundaries 
decreases as the angle of the boundary increases. The limiting size of the 
striations would then be given by the distance a vacancy can diffuse to 
cause this rearrangement before the mobility drops to too low a value. 
This model gives the correct order of size for the striations (Frank 1956 a) ; 
however, a more detailed picture for striation formation is necessary 
before the model can be substantiated. 


§ 5. THE INFLUENCE OF OTHER DISLOCATIONS 


In the preceding sections we have always assumed that, except for the 
dislocations formed by vacancy condensation, there are no other sinks for 
vacancies. However, as just discussed, there will generally be other dislo- 
cations grown in either by growth accidents or by the micro-segregation 
of impurities. Since these dislocations will eat up part of the excess 
vacancy concentration, the number of vacancies frozen-in at the interface 
will never be available for dislocation ring formation. No attempt will 
be made to discuss the annealing of vacancies in general and its compli- 
cated kinetics (Kimura ef al. 1959). We want only to estimate how many 
of the vacancies could possibly anneal out at the dislocations already 
present. For this purpose we assume again that a dislocation line forms 
a perfect sink for vacancies. This assumption may not be valid for metals 
in which the jog energy is high, i.e. in f.c.c. metals with low stacking fault 
energies. Specifically we wish to estimate how these grown-in, dislo- 
cations affect the nucleation of dislocation rings in order to answer the 
following questions : (i) at what relative undercooling, (AT /T'm)e, can dislo- 
cation rings grow as a function of the grown-in dislocation concentration 
and cooling rate, (ii) how many rings of a given area can be produced from 
the excess vacancies and what would be the resulting dislocation density 
from this source, (iii) if the rmgs should grow to form a network what 
would be the dislocation density? 
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It is possible to derive a differential equation for the actual vacancy 
concentration as a function of temperature, cooling rate, @, and grown-in 
dislocation concentration, pp, as parameters by considering the climb rate 
of these dislocations. However, if one wishes to study the solution as a 
function of these parameters, very tedious numerical integrations are 
necessary. Therefore, we calculate the excess vacancy concentration as 
a function of temperature by the following approximate method. 

We consider the continuous cooling as being approximated by a series 
of sudden quenches over small temperature intervals, A7;, with holding 
times of At; at each temperature. The cooling rate will then be given by 
Q=AT;/At;. In the time intervals At; the area 2;2 around a dislocation 
line drained of excess vacancies is of the order 7;2=D,At;. Let us now 
call the concentration of vacancies in thermal equilibrium after the ith 
quenching step c; (cj is identical with c, from eqn. 8) and the actual con- 
centration p;. Let us assume that the vacancies are uniformly distributed. 
If we have now py grown-in dislocation lines/em?2 they will drain an area 
Povi? per unit area in the time At; from the excess concentration (p;—ci). 
If we start with the equilibrium concentration cy) at temperature 7'y, then 
after the first quenching step the excess concentration, will be 


(P1— C1) =(Cgo—C))(1—poty?) . . . 2 ee), 
and for each successive step it follows that 
i Cie; — Ci) (1 ——pyti*)> ae 2 (20) 


By using the abbreviations €£,;=(1—%D,AT\/Tm), %=pol'm/Q and 
Aci =ci_, —¢; it can be easily shown with the aid of the recursion, relation 
(20) that 


Pr—en= ¥ ACm| TT &| - ee hi6 


i=™ 


after the cooling interval ATJ= > AT. 
i=1 

The equilibrium vacancy concentration, c,, and the vacancy excess, 
Pn—Cn; are plotted in fig. 5 as a function of AT7'/7'm with } as a parameter 
assuming the activation energies used for eqn. (17). From these curves 
we see that the final vacancy concentration, in the absence of disc 
formation, depends strongly on the parameter +. 

With the aid of fig. 5 we can now discuss the previously proposed 
questions : 

(i) The growth of vacancy discs will take place when the supersaturation 
PnlCn Yeaches a critical value Z given in table 3. The undercooling 
(AT /T'm)c is easily found in fig. 5 by noting that log p,, has to exceed logc,, 
by a fixed amount logZ. As it is readily seen (AT /Tm)c depends very 
weakly on the magnitude of Z except for very small values of ys. We take 
therefore as a typical value Z= 10, and have listed the resulting (ATT m)e 
in table 3 for the case 7’ = 1000°K, Q= 0-1°c/sec and various values of po. 
This cooling rate is representative for careful crystal growth so we see that 
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the presence of 105-10 grown-in dislocation lines per em? more than 

doubles the necessary undercooling (A7'/7'm)c for dislocation loop growth 

compared with a dislocation free crystal. 
(ii) The maximum number, N/cm* of dislocation loops of area, A, which 

can grow after the crystal has been cooled below the critical ratio (A7'/7'm)c 

is given b 

Aerie N=(p,—¢,)/bDA . a ee en 2) 


Fig. 5 


EOgmeG., 


©) , ae fo 


4 aS 6 af AS} a, 
AUP ns 
Plot used to determine the critical undercooling, (A7')7T m)e necessary to 
nucleate dislocation loops in the presence of pp grown-in dislocation 


lines/cm? as a function of the cooling rate, Q. (c,,=equilibriuam vacancy 
concentration, p,,—C,, = vacancy excess, f= pol’ m/Q.) 


where b is the Burgers vector. In order for these loops to act as Frank— 
Read sources during plastic deformation, the loops should be of the order 
of 1jin diameter. The number of rings for this radius r= 1 and also for 
r= 10, have been listed in table 3 and we see that an appreciable number 
of such sources could be produced during crystal growth over a wide 
range of grown-in dislocation densities. The resulting dislocation density 
p, is given by 

PR (=p DEE ah so ees ee ct, Se Cee 
and is also listed in table 3 (total dislocation density, p’=py+p). Loops 
of radius r = 1 p cannot be detected by etch-pit techniques, and for r= 10 p, 
p must be about 10°cm~ to appear uniformly distributed by etch-pit 
observation. Thus, if py is 104-108cm~? a uniformly distributed array of 


dislocation etch-pits could not be due te dislocation loops formed by 
vacancy condensation. 
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(iii) The loops can grow to a network, if their radius r grows to a size 
comparable with their spacing, i.e. if N=1/r?. Thus the dislocation 
density, pnet, in a plane cut by the network is given by 


1 
pnet® 5 © (Pn — en)» /b*. eee han Cee) 


For the particular case of 7=1000°K and Q=0-1°c/sec, pnet is also 
tabulated in table 3 for various pp. 

From these considerations it is not possible to tell whether the 
dislocations found in crystals are due to vacancy condensation or due to 
other mechanisms. The only conclusion that can be reached is that if in 
slowly grown crystals, a dislocation density pys 10°em™ of grown-in 
dislocations exists, then there are not enough vacancies left to form an 
appreciable network. However, with fast quenching rates from high 
temperatures the crystal can contain a large number of dislocation loops 
formed by vacancy collapse even when the crystal contains appreciable 
numbers of grown-in dislocations. 

Finally, we want to discuss very briefly some experiments on highly 
perfect Si crystals. Dash (1958) has postulated that dislocation climb in 
the solid behind the advancing solid—liquid interface allows the dislo- 
cations to move to the outer surface of his crystals; thus leaving them 
dislocation free. However, as we have already shown, the climb rate is 
extremely slow relative to the rate of advance of the interface at small 
distances from the interface. Thus, this mechanism does not seem very 
likely to have produced the dislocation free crystals. It is much more 
likely that the dislocations were removed from his crystals due to the 
shape of the solid—liquid interface during a portion of the crystal growth 
(interface convex to liquid). 


APPENDIX I 


ATTRACTION OF VACANCIES TO DISLOCATIONS DUE TO 
Evastic INTERACTION 


The interaction potential of vacancies with dislocations is (Friedel 1956 b) 
GeW fo 8) aan oe 


where W is the interaction constant, ) the interatomic distance and r the 
distance from the dislocation line. Thus, the drift velocity v of vacancies 
towards the dislocation due to the interaction potential is 
fe Row, 2WDt 

kp or aaa kT rs . . . - . . . (A 2) 
The time f) it takes the vacancies within a radius rg to anneal out at the 
dislocation, due to the drift velocity in the energy gradient, is then 
easily shown to be 

BD 


i= SWDb* er ee re Oe ne 8 (A3) 
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From eqn. (14) we use 7)7=b/7AC as the radius of the region which can be 
drained by the dislocation before it moves out of the drained region. 

By the same argument as used in the text we obtain then as the cli mbing 
velocity 

ae Loess WDAC® gids sopra tase CA) 
bo bkT 
This velocity vp’ would be obtained if vacancies would move towards the 
dislocation exclusively due to the action of the elastic attraction potential. 
With typical values W z 0-5ev, 7T'~ 1000°K and AC Zz 10-3 this velocity is 
always several orders of magnitude smaller than the steady-state velocity 
obtained from eqn. (15) due to chemical diffusion. 

That we can neglect the elastic attraction of vacancies in the climbing 
process can also be seen by the following argument: The average chemical 
potential gradient during the steady state climbing at temperature 7’ is of 
the order of :/7) where yu is the chemical potential of vacancies quenched- 
in at the temperature 7’. If the sample is cooled by a temperature 
interval A7’ we obtain 


pe —wO;ATAC 

ro a bT m 
(U;=energy of vacancy formation, AC = excess concentration over thermal 
equilibrium value, b=atomic distance). The elastic energy gradient is 
however from eqn. (A 1) 


(A5) 


if ¢ 2 
AL as ile (A 6) 


or ry? 
By equating (A 5) and (A 6) we obtain the radius 7, at which both gradients 
have the same value. It can be shown that in all cases r, is less than a 
few percent of 7). Hence the attraction potential can be neglected over 
the greatest part of the diffusion path. In typical quenching experiments 
r, is of the order of 100 or less. The presence of the elastic attraction can 
then be accounted for by replacing the dislocation by a sink of effective 
radius 7, + 100. 


APPENDIX II 
DISLOCATION CLIMB RATE 


We assume a straight dislocation line lying in the z-axis of a cartesian _ 
system (x,y,z). It acts as a sink for vacancies and may move with 
constant velocity v, in the x direction. We also assume that a steady 
state is reached in the coordinate system (g, y,z) moving with the 
dislocation. vi 

If r is the distance from the dislocation line, Dy is the diffusion 
coefficient for vacancies and cm the concentration of vacancies at infinity, 
the concentration of vacancies at the point (r, €) is given by Rosenthal 


(1946) to be 
0) Vol 
em—e=Aexp(— 35) Ko(sh), Mare. ae (AT) 


where K, is the modified Bessel function of second kind, zero order. 
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The constant A will be determined as follows: As discussed in Appendix 

I, we use rs= 100 as sink radius for the dislocation. The velocity of the 

dislocation, line is proportional to the flux of vacancies over a cylinder of 
radius rs around the dislocation line, or 


= 2Q7rDyrs (5) = a VoD. . . . . . . (A 8) 


& 


If rs< Dy/v, we obtain then 


_ ad) ee 
ior i exp( re) (3). ) 


Now we have the additional condition that at r=r, the concentration c 
has its thermal equilibrium value cy). From eqn. (A9) we obtain, then, 
an, equation to determine the velocity vy as a function of ¢m and cy which 
leads, for rs << Dy/vp to 


Qqr Dy(€m— C9) 


= —____—_.. S ett ee sak) 
In (4D,/9v975) b 


Vo 


If we compare the exact result (A 10) with the velocity given by eqn. (15) 
we find that eqn. (15) must be multiplied by a factor f= 2/In(4D,y/Svprs), 
where $ is Euler’s constant (9=1-78). Since all variables in f give only a 
logarithmic dependence we may take average values for Dy and vp» 
(Dy x 10-> cm?/sec, rg ¥ 3 x 10-7, v9 = 10-3 cm/sec) and obtain f= 0-2 to 0-1. 
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ABSTRACT 


The object of this paper is to estimate, semi-empirically, the Fermi energy of 
the free electrons in a metal and to indicate why a simple correlation may be 
expected between this quantity and the electronic work function. Such a 
straightforward correlation is found to hold approximately for most metals in 
Groups I to IV of the periodic table. 


§ 1. INTRODUCTION 


Tue theoretical considerations of Wigner and Seitz (1933, 1934) on the 
cohesive energies of the alkali metals suggest that the energies of the 
electrons in the metal may be described to a fair approximation in the 
following manner. The valence electrons move under the influence of the 
periodic electrostatic field of the ion cores and also under their mutual 
interactions. In accordance with Fermi—Dirac statistics the electrons 
occupy a band of energy levels, based on a ground state energy, Hy, up 
to a maximum energy, (H,+ Co), at absolute zero. EH, is the lowest energy, 
potential plus minimum kinetic energy, which a vclence electron can 
have in the field of a single positive ion core within the metal lattice, with 
the appropriate boundary conditions. ( is the additional kinetic energy of 
the most energetic electrons, since, by the requirement of Pauli’s exclusion 
principle, all electrons cannot occupy the same quantum state. f, may 
be assumed to have the same value as for free electrons moving in a bounded 
region, of constant potential V5, outside which V is zero and within which 
the electron density is the same as in the metal under consideration. 

In this wave mechanical treatment the mutual electron interactions 
give a further large positive contribution to the total energy not already 
accounted for in Hy. This is offset by a negative contribution, arising from 
the exchange interactions between electrons with like spin and different 
wave numbers. Wigner and Seitz (1933, 1934) calculated the positive 
term on the assumption of a uniform distribution of negative charge and 
this was corrected to allow for the correlated motion of the electrons 
(Wigner 1934). The effect of this correlation and also of the exchange 
energy is to tend to keep the electrons away from each other and thus to 
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reduce the mutual electrostatic interactions of a uniform electron distri- 
bution. The calculations indicate that together, for free electrons, these 
three contributions almost cancel one another, leaving only a small net 
interaction energy, Q;, of the order of 10% of the average energy of an 
electron in the metal relative to an electron at rest in free space as the 
energy zero. 

The mean energy of an electron in the metal is thus approximately 

H= Hy + 2f, 
where the factor 3 is the averaging factor relating the mean additional 
kinetic energy of a free electron distribution to the maximum additional 
kinetic energy, Co. 

In the absence of any double layer at the metal surface least work is 
required to remove an electron from inside the metal to rest at infinity 
if an electron of maximum energy is ejected. This will require an amount 
of work of the order — (H+ ¢,) at absolute zero if there is no relaxation of 
the electronic charge distribution remaining within the metal. There will 
in fact be some relaxation of the electrons within the metal when, one is 
ejected and therefore the energy required to remove one electron will be 
less than —(H,) +9) by the amount of this relaxation energy (— AQ,) say. 

Consider a ring of metal which is cut to create fresh surface. At the 
instant of separation, of the two surfaces the electrons associated with the 
new surface ions have their potential energy increased (less negative) to 
approximately half the numerical value which was obtained before 
separation. ‘This is because the electrons instantaneously associated with 
the surface ions are no longer completely surrounded by positive ions. 
The kinetic energy of the electrons in these surface states will, by the 
virial theorem, be reduced by a factor of approximately one half also, so 
that the numerical magnitude of the total mean energy of an electron in 
a surface state will be of the order of one half that of the mean energy of 
an, electron in the interior of the metal. Similarly, the magnitudes of the 
maximum electron energies in the interior (H,)+¢)) and in the surface 
(EZ, +), will be in the approximate ratio of two to one. 

Electrons will therefore move from the surface atomic cells toward: 
the interior of the metal setting up a double layer of charge, with the 
positive layer on the outside. Equilibrium will be established when the 
work done to move an, electron against the double layer field from the 
surface to the interior of the metal is equal and opposite to the decrease 
in energy (Ey + %o) —(Zy+%)s- The work, x, required to remove an, electron 
from anywhere in the system, interior or surface, to infinity is now the 
same. Thus to remove an electron from the interior requires an amount of 
work : 

y= —(Eyoth)—4ane’a . . . . - ss (I) 
where 47ne2x is the work gained in passing through the double layer of 
strength n electrons/sq. cm and thickness 2. But 


Agneta + (Ey +6) — (Hot So)s=0 + » © + + (2) 


P.M. E 
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Therefore y in (1) also equals (Z)+ 9), which is the work required to 
remove an, electron from the surface layer. 

The work required to remove an electron at constant temperature and 
pressure from the system is a measure of the change in the free energy of 
the system per electron removed, which is the negative of the thermo- 
dynamic potential of the electrons. At equilibrium therefore the thermo- 
dynamic potential of the electrons in the metal is everywhere the same, 
as indeed it must be. For other and further discussions of electronic work 
functions reference may be made to the review article by Herring and 
Nichols (1949) and also the book by Seitz (1940). 

The foregoing argument suggests that the work function y should be 
of the order of —4(H,+¢,), if the mutual interaction terms within the 
electron distribution of the metal can be neglected. 

In the next section semi-empirical calculations of the Fermi energy 
C=(£,+%,), are made using as many experimental facts as possible, and 
assuming that the resultant electron—electron interaction energy can in 
fact be neglected within the metal, but not between valence electrons in 
polyvalent atoms. The arguments used in computing ¢ for elements with 
more than one valence electron have been set down, before by the author 
(Varley 1954) but are reproduced here for completeness. 


§ 2. EXPRESSIONS FOR THE FERMI ENERGY OF ELECTRONS IN METALS 
2.1. Monovalent Metals 


The average energy of an electron in the metal, relative to the energy 
zero of an electron at rest in free space, is taken to be (H,+3¢,). The 
electrons within the metal can be brought to this zero energy state by 
vaporizing the solid and then ionizing the atoms. It follows that 


Ey+30,= —(A+J,) . e . . . . . (3) 
where A is the latent heat of vaporization per atom at absolute zero and 


I, is the first ionization energy of the free atom. 
Thus #y)= —(A+1,+ 26), 


t h2 3n\ 2/8 
Tairen (*) ; 
where is the number of free electrons per unit volume, h is Planck’s 
constant and m is the electronic mass. The Fermi energy for monovalent 
metals is thus: 
h? (3n\28 
C= (By+)=—| A+L,—- soa(+) i Se 


20m \ a 


2.2. Polyvalent Metals of Groups II, III and IV 


In a polyvalent metal the exchange interactions between electrons of 
like spin approximately cancel out the repulsive, correlated Coulomb 
interactions. This effect of the exchange forces is lost when the metal is 
vaporized into free atoms and in the latter the repulsive electrostatic 
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forces predominate. The exchange forces tend to keep electrons of like 
spin apart so that most of the exchange energy may be regarded as arising 
between electrons which are in different atomic cells. Therefore the 
energy per atom required to vaporize a polyvalent metal into free atoms 
may be considered to result from (a) changing the energy of the Z valency 
electrons associated with the ion core of each atom and (bd) overcoming 
the exchange forces which produce an effective attraction between adjacent 
atomic cells. 

The mean energy of an electron in the metal is (H)+ 2); its energy 
in the free atom is 


Il @ 
=5 21 
Zi 
where J; is the ionization energy of the ith valency electron and the sum 


is over the Z valency electrons. The change in energy on evaporation is 
thus : 


is 2 
a7 2 Fir (Bo+ $60). 


The exchange interaction energy between electrons of like spin in different 
atomic cells of the metal is approximately the difference between, the 
total interaction energy per electron in the metal and that in the free 
atom. The total mutual interaction energy per electron in the metal is 
taken to be zero. That in the free atom is 


LZ Z—2 
{lem 3b lp 

where J, is the ionization potential of the last valence electron and Jp is the 
energy to promote an electron from the s to the p state in the ionization 
spectrum M+. The last term only occurs in elements of valency greater 
than 2. 

The energy of vaporization per electron is A/Z, where A is the experi- 
mental heat of vaporization per atom, and this is therefore equal to 


i ema BE 
Fila GU (Bot 8). + + 


= 


It follows that 
Z—2 A h2 (=)" (6) 


2 Z 
b= (Hot b)=— 7% litte Ft 7+ 30m 


for an element with Z free electrons per atom. 


7 


§ 3. NuMERICAL COMPUTATIONS 


Using data from standard references and published papers (H andbook 
of Chemistry and Physics 1957-8, Kubaschewski and Evans 1951 ) on ioniza-— 
tion spectra and heats of vaporization, calculations of the Fermi energy, ¢, 
have been made for several metals in the first four groups of the periodic 
table. These calculated values, together with experimental values of the 
work function, y (Handbook of Chemistry and Physics 1957-8) are set out in 


E2 
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the table together with the ratio (—x/¢) in the last column. It will be 
seen that in most cases this ratio is not too different from the value 0-5 
suggested in the text. 


Element —¢ x =e 
Li 5-08 2-49 0:49 
Na 5-01 2:28 0:46 
1K 4-45 2-24 0:50 
Rb 4-34 2-09 0-48 
Cs 4-06 1-81 0:45 
Cu 8°38 4-5 0-54 
Ag 8:33 4-73 0-57 
Au 10:89 4-82 0-44 
Be 6:88 3-9 0-57 
Mg 5:50 3-68 0-67 
Sr 5:78 2°74 0-47 
Ba Dehe, 2°5 0-44 
Zn. 6-22 4-0 0-64 
Cd. 6-54 4-07 0:62 
Hg 7:85 4-5 0:57 
Al 5-64 4-0 0-71 
Ga 6-75 4-] 0-61 
aN 7-96 3-68 0-46 
Ge 7:99 4-5 0-56 
Sn etd 4:38 0-56 
Pb 8-13 4-() 0-49 


§ 4. ConcLUDING REMARKS 


An exact ratio of 0-5 for (—y/¢) is, of course, not to be expected for 
several reasons. There are many approximations in the caleulations for 
the Fermi energy. For instance, it has been assumed that the energy to 
ionize an electron in a free monovalent atom gives a measure of the energy 
of that electron alone, whereas the electrons in the ion core will relax 
and their energies will be reduced when the valence electron leaves the 
atom. 

Again, using the free electron value of ¢) is only approximately true, 
although soft x-ray spectra indicate that the true energy band widths 
are not very different from the free electron widths in several cases. How- 
ever, differences are to be expected, especially in the case of the polyvalent 
elements in which the details of the energy bands are influenced by the 
periodic lattice potential (Brillouin zones) at the top of the Fermi distri- 
bution. Nevertheless, an approximate correlation between Fermi energy 
and work function seems to be indicated on the basis of these calculations 


and may possibly serve as a rough guide for other elements, such as the 
transition elements, in the periodic table. 
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ABSTRACT 


Refinements of the curve of density of states n(#) in lead, as deduced from 
the de Haas—van Alphen effect, are inferred from existing superconductivity 
and transport data on pure lead and on Pb—Tl and Pb—Bi alloys when these 
data are interpreted in terms of a simple ‘inverse transition metal ’ 
model. (£) is found to have a maximum just below the Fermi level. It is 
suggested that the maximum might be accounted for by supposing the 
multiply-connected portion of the Fermi surface to split up into separate 
portions as the electron concentration is reduced; however, this particular 
interpretation does not seem to be compatible with a model of rigid bands. 


A HIGHLY schematic curve of density of electronic states n(£) versus 
energy H for lead has been inferred recently from an interpretation of the 
de Haas—van Alphen effect (Gold 1958; hereinafter referred to as 1) and 
is shown in fig. 1. Also shown are the contributions from the various 
Brillouin zones; the curves have been, drawn in a more or less parabolic 
fashion from the computed values of n(H,) at the Fermi level such that 
their areas give the number of electronic states within the various parts 
of the Fermi surface. In order to obtain a more detailed n(#) curve for 
pure lead we shall examine first of all the effects of varying the electron 
concentration as inferred from existing data on, alloys of lead with thallium 
and bismuth. Since thallium and bismuth occur on either side of lead in 
the same row of the periodic table, alloying with these particular metals 
should produce the least distortions in the band structure of pure lead. 
The variation of the density of states n(H,) at the Fermi level with 
decreasing electron concentration can, be obtained from the superconducting 
critical-field curves of Pb—T] alloys; no suitable data are available for the 
Pb-Bi system. n(H ) is related to the electronic specific heat constant. 

yei by 
n( Eo) = 1-59 x 10-19(a3/4)(3/277?k?) yer per ev per atom (1) 


(cf. Mott and Jones 1936, p. 179), where a is the lattice constant in em 
and ye) is given by 
Ye = 1/27(H,/T',)? erg om-3 deg-*, ae eg 


+ Communicated by the Author. 
{ National Research Laboratories Postdoctorate Fellow. 
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assuming the critical-field curves to be parabolic (Shoenberg 1952, Ch. IIT). 
The mean, zero-field transition temperatures T, have been interpolated 
from the results of Meissner et al. (1932) and the critical fields H, at the 
absolute zero have been extrapolated from the data of Schubnikow et al. 
(1936), taking the critical field to be the value at which magnetic flux 
begins to enter the specimen rapidly as the field is increased. The 
dependence of a on Tl concentration has been given by Tang and Pauling 
(1952) and the resulting variation of n(#)) deduced from (1) and (2) is 
shown in fig. 2(a). Also shown is the value of n(Z,) obtained by Decker 
et al. (1958) from the actual non-parabolic critical-field curve of pure lead ; 
however, the data of Schubnikow et al. are insufficient to show up any 
significant departures from parabolicity for the alloys. 


Fig. 1 
1.2 
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Schematic curve of density of states near the Fermi level as inferred from the | 
de Haas—van Alphen effect (I). The energy is measured relative to the 
Fermi level. 


While this interpretation of the superconductivity data is likely to be 
too naive on account of the lack of a Meissner effect in the alloys, the 
curve of fig. 2(a) nevertheless suggests that n(H)) drops rapidly fom its 
value for pure lead, reaching some fairly constant value at about 15 HxAlily 
It is worth noting that the yei values obtained by Stout and Guttman (1952) 
from the critical-field data on alloys of tin with up to 20% thallium show 
no such rapid variation, thus suggesting that the decrease found in the 
Pb-T] system is indeed an electronic structure effect. 
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We can obtain further information about the variation of n(£y) on 
alloying from the temperature dependence of the resistivity, 0p/oT, and 
from the absolute thermoelectric power S. In interpreting such data it 
is convenient to think of lead as a kind of ‘inverse transition metal’ since 


(a) 
04 n(Eo) 
per eV per atom 


9) 
_— oe- 
5 10 f° 
(0) 2 oO 10% x (3) 
o°c 


8 aT, 

A pp 2 cm. deg-' 
6 
4 


(c) 


n(E,) 


(4) ~0.6. 
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50 40 30 20 10 ie) ite) 20 


awa cout % Bi 


Variation with Tl and Bi concentration of (a) n(H) from critical-field data, 

(0) p/@7' at 0°c, (c) S at 0°o and (d) n(£y) asinferred from (a), (6) and (c). 

Curve (a) is from Meissner et al. (1932) and Schubnikow et al. (1936); the 

Pb-Tl portion of curve (b) is from Guertler and Schulze (1923) and 

Tammann and Riidiger (1930) and is estimated to be reliable to about 

2%. @, Decker etal. (1958), assuming parabolic H, vs. 7 curve; O, Decker 

et al., from actual non-parabolic H, vs. T' curve; @,. Vv, Thomas and 

Evans (1933); A, Rolla (1915); (9, Christian et al. (1958). The a-phase 

boundaries lie at 80% Tl and 18% Bi (cf. Hansen 1958), with ordering 
taking place at PbTls (ie. 75% Tl; Tang and Pauling 1952). 
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we see from fig. 1 that the rapidly varying contribution from the electrons 
to n(#) is much greater than the slowly varying contribution from the 
holes. On the basis of such a model we should expect the major part of the 
conductivity to be due to the holes (cf. Mott and Jones 1936, -p. 267); 
this has been confirmed experimentally by Borovik (1954) who has shown, 
on. the basis of a two-band analysis of the galvanomagnetic effects, that the 
holes in fact account for between 85%, and 90°% of the total conductivity 
in pure lead. Thus, assuming the validity of the Mathiessen rule, we may 
write for the temperature dependence of resistivity (cf. Mott and J ones 
1936, p. 301) 

Op Opp  O 

a? oT OT a 
where Wy is the number of holes per cc; we have taken 1 [tr <Tn( Eo) 
and have assumed that the effective mass of the holes is not appreciably 
altered on alloying (mn*(Pb)=1-0m,) from I). Turning now to the 
thermoelectric power, we have rather generally 


2227) C 
Ree eet 2108 LI) | oa sh vot ae ann (4) 

3e oH 0 : 
for our model this reduces to (cf. Mott and Jones 1936, p. 313) 


joe em es) 
0 


(Ho) 


(mn*/. f eT OC j 
NW h 


3e n(H,)\ OH 
1 /dn(E) 
=2 2°45 7 —___ [| = - oe, 
- Ey ( an ),e eae a 


the energy being measured in ev. 

The variation of dp/dT' at 0°c is given in fig. 2(b); the Pb—Bi results 
are from Thomas and Evans (1933) while the Pb—T!1 curve has been drawn 
from values of the temperature coefficient of resistance given, by Guertler 
and Schulze (1923) and of the absolute resistivity given by Tammann and 
Riidiger (1930) (the absolute resistivity measurements of Guertler and 
Schulze are fraught with systematic errors). As we move from left to right 
in fig. 2(b) the number of holes in the second zone will be progressively 
reduced so that we should expect from (3) a smooth increase in dp/d7’, 
modulated by n(E,). We therefore infer that n(H#)) begins to increase 
from a roughly constant value near A, reaches a maximum near B and 
then falls through its value for pure lead, this decrease being offset near C 
by the effect of diminishing 1); the flattening-off of the curve at D is 
probably due to an increasing contribution to the conductivity from the 
electrons so that (3) is no longer valid. Thus the behaviour of dp/d7 for 
the Pb—TI series confirms qualitatively the scheme of fig. 2(a) but it is 
difficult to draw any definite conclusions from the Pb—Bi results. 

In fig. 2(c) we show the variation of S at 0°c from Rolla (1915) and 
Thomas and Evans (1933), the relative values given by these authors 
having been converted into absolute values using the absolute value of S 
for pure lead given by Christian et al. (1958). Although the variation of 
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S with 7 for pure lead is approximately linear above about 100°K ( Christian 
et al.), it is not exactly in accordance with (5) since the high-temperature 
results do not extrapolate to S=0 at 7=0; this matter will be discussed 
later but for the present we shall take (5) to be approximately valid for 
our discussion of the alloys. While the effect'of impurity scattering on 
S is not well understood we shall nevertheless suppose it to be small by 
comparison with the rapid rise in S with increasing Tl concentration. 
Assuming then a quasi rigid-band model for the alloys it seems reasonable 
from (5) to associate the large positive S at 6% Tl with a sharply rising 
n(E,) and the value S=0 at about 2-5% Tl with a maximum in n(£p); 
both inferences are in qualitative agreement with fig. 2(a). The negative 
S near pure lead implies a falling n(H)) while the subsequent slow rise 
in S with increasing Bi concentration suggests that n(#)) begins to 
flatten-off somewhat; in view of the inapplicability of our model when the 
electronic contribution to the conductivity becomes appreciable, the rising 
n(H,) beyond 11% Bi as suggested by (4) must be viewed with some 
scepticism. 

We are now in a position to give, in fig. 2(d), an approximate curve of 
the variation of n(#,)) with Tl and Bi concentration. For pure lead we 
have taken the value from Decker et al. as being the most reliable. Although 
the curve of fig. 2(a) suggests a total drop in n(£,) of roughly 60% with 
the addition of about 15°% Tl we shall take this drop to be nearer 30% 
as being more in keeping with the results of fig. 2(b). Assuming then that 
any departures from the rigid-band model are not too large, so that changes 
in n(£,) on, alloying can be related to the n(£) curve for pure lead solely 
by a shift in the Fermi level, we give in fig. 3 a more realistic density of 
states curve for pure lead than that given in fig. 1, changes in area under 
the n(£) curve being related to the corresponding changes in electron 
concentration in fig. 2(d). The peak in the n(£) curve is then found to 
lie at about 0-02ev below the Fermi level. The contributions from the 
holes and overspill electrons are also indicated; as is required for a metal 
of even valency, these have been drawn in with equal area, this corres- 
ponding to a weighted mean of the number per atom of holes and electrons 
given in fig. 1, viz. 0-41+0-06 per atom. 

We shall now give a discussion of the temperature dependence of the 
thermoelectric power of pure lead (Christian et al. 1958) in relation to the 
n(H) curve. The low-temperature values of S below about 6°xK lie on a 
straight line with S=0 at 7'=0 (results of Steele; see fig. 3 of Christian 
et al.) and we have used these results and (5) to infer the slope of the n(Z) 
curve at the Fermi level (broken line in fig. 3); in this connection we note 
that the variation of S with addition of small amounts of impurities is in 
no sense ‘anomalous’ (Jan et al. 1958). As we have already mentioned, S 
also varies linearly with 7’ above 100°K, but the variation is much slower 
than that below 6°K and, moreover, the high-temperature results do not 
extrapolate to S=0 at 7'=0 in accordance with (5). These features may 
be accounted for qualitatively in terms of the rapid change in On(E)/0F . 
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near the Fermi level from positive to negative values within an energy 
range comparable with kT. We can make a crude estimate of this effect 
by considering the n(#) curve near the peak to consist of two intersecting 
straight lines and by supposing that we can replace (On(£)/0F), in (5) by 
a straightforward average of dn(H)/dE over an energy range of 2k7' on 
either side of the Fermi level. This results in the modification of (5) by 
the factor ((C'+ 1)«/2k7’—C+1)/2, where C is the magnitude of the ratio 
of dn(H#)/0H below the peak to that above the peak and « is the energy 
between the peak and the Fermi level (2k7'>.). The modified linear 
expression for S can be fitted to the high-temperature results by choosing 
C~0-8 and «~0-0lev. This crude argument thus suggests that the slopes 
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Proposed curve of density of states. The energy is measured relative to the 
Fermi level; ——-, slope at the Fermi level from absolute thermo- 


electric power below 6°K ; O, Decker e¢ al. (1958). 


of the n(Z) curve on either side of the peak are comparable in magnitude 
and the value of ¢ turns out to be comparable with that suggested by our - 
rigid-band interpretation of the alloy resultsf. 

It is tempting, at first sight, to account for the peak in the n(Z) curve 
at roughly 0-02 ev below the Fermi level in terms of the multiply -connected 
Fermi surface in the third Brillouin zone, a portion of which is illustrated 
diagrammatically in fig. 4(a); a drawing showing the three-dimensional 
character of this surface is given in I. When the electron concentration 
has been, sufficiently reduced (e.g. by alloying with Tl), this surface will 
break up into separate portions centred about either the zone corners W 


+ We should also expect the electronic contribution to the specific heat 
for I> ~100°K to be much less than that deduced from low-temperature 


measurements. 
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(fig. 4 (b)) or the mid-points K of the lines joining adjacent corners (fig. 4 (c)), 
depending on, whether the energy is lowest at W or at K respectively. 
When, the separate portions are nearly in contact, grad, # along zone 
lines will be particularly low, resulting in a peak in n(£) when contact 
has just taken place. However, if this is the correct explanation, it is 
indeed surprising that the peak should lie so close to the Fermi level. 
The ‘limbs’ of the multiply-connected surface are, in fact, quite thick 


Fig. 4 


Diagrammatic illustration of the breaking-up of the multiply-connected 
third-zone Fermi surface. 


so that one would expect the breaking-up of the surface to occur at an 
energy comparable with that associated with complete disappearance of 
the surface, i.e. at roughly lev below the Fermi level; thus the above 
interpretation of the peak does not seem to be compatible with a model of 
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rigid bands. In order to assess the extent of shifts in the band structure, 
low-temperature heat-capacity measurements on isoelectronic Pb—TI-Bi 
alloys would be very desirable (cf. Guthrie 1959). 

In conclusion, we wish to point out that there is some doubt as to the 
existence of the cushion-shaped pieces of Fermi surface in the fourth zone. 
From the detailed reasoning given in I we can state that either there are 
electrons in this zone giving a reasonably large contribution, with the 
lowest energies in zones three and four being nearly equal (fis. 1) ,20P 
the fourth zone is completely empty; the evidence from the de Haas—van 
Alphen effect alone is insufficient to distinguish unambiguously between 
these two possibilities. However, the values of the electronic specific 
heat and of the total area of the Fermi surface as deduced from the de 
Haas—van Alphen data can be brought into considerably better agreement 
with the results of other experiments by taking the fourth zone to be 
empty; this supposition also makes the numbers of electrons and holes 
more nearly equal. Since in drawing fig. 3 the data in fig. 1 for the holes 
have been more heavily weighted than that for the electrons, it makes little 
difference to the areas under the curves whether a contribution from the 
fourth zone is included or not. 
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ABSTRACT 


High resolution electron microscopy has been used to examine thin films of 
natural rubber and other polymers. Dispersions of various materials (in 
particular neoprene and chlorinated natural rubber) in natural rubber have 
been studied. Structures which are believed to be the individual curled-up 
molecules of the dispersed polymer have been photographed. Stretched 
rubber films have also been examined and some information about the mode 
of deformation obtained. 


§ 1. LyTRODUCTION 


Frims of various plastic materials have been commonly employed in 
electron, microscopy—for the most part as supports for other specimens 
and as replicas of solid surfaces from which they have been stripped. 
They have not, however, been used very much as specimens in their 
own right to provide information about the properties of the plastic 
material itself (a notable exception is the work of Agar et al. 1959, on 
polyethylene). In the present paper an account will be given of an electron 
microscope study of rubber films from which information about the rubber 
has been deduced. In particular, structures which are thought to be 
individual rubber molecules have been found. 


§ 2. MATERIALS AND TECHNIQUE 
2.1. Materials 


Films of four materials were examined. These were natural rubber and 
three related polymers—neoprene, guttapercha, and Alloprene. Natural 
rubber is mainly cis-polyisoprene (—CH,—C . CH, = CH—CH,—),,, neo- 
prene is polychloroprene (—CH,—C.Cl=CH—CH,—),, guttapercha is 
mainly trans-polyisoprene, and Alloprene is a chlorinated natural rubber 
made by Imperial Chemical Industries Limited (approximate composition 
CioH,3Cl,). A purified natural rubber was also used. This was rubber 
from which the small amounts of proteins, ete. which commonly occur in 
the natural product had been removed. None of these materials had been 


y+ Communicated by Dr. F. P. Bowden, F.R.S. 
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vulcanized. Thus, there was little cross-linking of the molecular chains 
and the polymers were soluble in suitable organic solvents. All the samples 
were obtained from the British Rubber Producers’ Research Association. 


2.2. Specimen Preparation 

A solution in ‘Analar’ benzene of the material to be examined was 
made. A drop of this was allowed to fall on a clean, water surface where 
it spread and, after evaporation of the benzene, left a floating film. Part 
of this was then picked up on a specimen grid and subsequently examined 
in the electron microscope—Siemens Elmiskop I. The film thickness was 
found to be critical. The concentration of the solution was therefore 
adjusted to give films of the correct thickness judged by the interference 
colour—low order white with veining of yellow. Such films (of the order 
of 1000 A thick) were over most of their area too thick for high resolution 
electron microscopy, but they did show thin patches where the local 
thickness was probably less than 1004. These were the regions studied. 
Films of lower mean thickness could not be used since they formed networks 
of thin threads rather than continuous sheets. 

Stretched rubber films were also examined. The technique of stretching 
them was as follows. The film while floating on the water surface was 
touched with two pieces of filter paper. Unvulcanized rubber being sticky, 
it adhered to them. The pieces of filter paper were then drawn, apart and 
parts of the stretched film in between picked up on specimen grids. Once 
on a grid the film remained stretched since it adhered to the bars of the 
grid. 


§ 3. RESULTS 
3.1. Changes Produced in Specimens by the Electron Beam 


Polymers suffer changes when irradiated with electrons (see, for example, 
Little 1954, Charlesby and Pinner 1959). Consequently, it was thought 
advisable to investigate the effect on a specimen of its irradiation during 
examination in the electron microscope. No visible changes occurred in 
rubber films as the electron, beam intensity was increased from a low value, 
but a rubber film after examination in the electron microscope was found 
to be insoluble in benzene. This change was studied further by placing 
glass plates coated with a film of rubber in the camera of the electron 
microscope instead of the normal photographic plates, and exposing them 
to the much less intense beam of this part of the instrument. It was found 
that exposure for about 2 min to a fairly intense beam rendered the 
rubber insoluble in benzene. This suggested that, even if examined at 
low magnification (say 5000 x) and low intensity, a rubber film in the 
normal specimen position would be rendered insoluble within a few 
milliseconds. Consequently, if a rubber specimen is examined in the 
electron microscope one cannot avoid altering it by the action of the 


electron beam. 
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However, sharp spotty diffraction patterns were obtained from crystalline 
polymers (polyethylene, guttapercha) provided that a reasonably low 
intensity beam was used. Thus the changes in these polymers caused by 
a weak electron beam were not too drastic. It seems likely that the first 
effect of the electron beam is a cross-linking of the molecular chains, 
followed after much greater irradiation by a complete disruption of the 
structure. The rubber specimens if examined at low intensity probably 
remained as rubber (albeit vulcanized rubber), but even if an intense beam 
was used the changes in structure produced no visible alteration of the 
specimen. 

3.2. Appearance of Rubber Films 

An electron micrograph of a purified natural rubber specimen is 
reproduced in fig. lf. Such micrographs show little structure apart from 
gross features such as holes, wrinkles, thick and thin patches, etc. There 
is a granular structure on a scale of about 20 A visible at high magnification, 
but this may well be due to electron-induced contamination. Films of 
ordinary natural rubber, neoprene, Alloprene, and guttapercha have a 
similar appearance. 

Filled rubbers have also been examined and the dispersion of filler 
particles photographed (fig. 2). These are important since carbon black 
and other materials are commonly incorporated in rubber to improve its 
mechanical properties. A micrograph of a stretched filled rubber is shown 
in fig. 3. The thin patches (which occur sparsely in unstretched rubber) 
are more extensive and seenrto have developed at the expense of the thicker 
regions. This would be expected for unvuleanized rubber will flow 
plastically but the sharpness of the boundary between thick and thin 
regions is perhaps surprising. 


3.3. Mixed Films 


The most interesting results were obtained from films made from a 
solution, containing a mixture of polymers. Mixtures of purified natural 
rubber and neoprene, and purified natural rubber and Alloprene were used. 
The neoprene and Alloprene which contain a high percentage of chlorine 
might be expected to scatter electrons more strongly than the natural 
rubber which apart from impurities is entirely made up of carbon and 
hydrogen. Thus the neoprene or Alloprene should appear darker than the 
natural rubber in a bright field image. Films of natural rubber containing 
a few per cent of neoprene show numerous dark patches (fig. 4): these 
presumably are the neoprene. Furthermore films of neoprene containing 
a few per cent of natural rubber also have the expected appearance of 
light patches on a dark background (fig. 5), but it is difficult to distinguish 
these from the occasional thin patches which these films show. 

A more detailed study was made of the Alloprene/natural rubber system. 
Figure 6 shows the appearance of a film containing equal parts by weight 


+ With the exception of fig. 10, all figures are as plates. 
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of the two materials. The contrast system is difficult to interpret, but the 
film is obviously not homogeneous. One would not expect it to be, for 
chemically different polymers are, in the absence of a solvent, generally 
immiscible or only slightly miscible (Dobry and Boyer-Kawenoki 1947, 
Scott 1949). 

Figure 7 was obtained from a film containing 20% of Alloprene. The 
dark spots presumably represent globules of the Alloprene. The profiles 
of these globules could be seen in regions where the edge of a film had 
curled over. They were thicker than the surrounding film and appeared 
to be flattened spheres. Thus at least part of the contrast is here due to 
the greater thickness of these globules. Again one would not expect 
miscibility of the polymers in these proportions. Similar pictures were 
obtained from films containing 10° of Alloprene. A micrograph of a 
stretched film having this composition is shown in fig. 8. It can be seen 
that the film has torn at the Alloprene globules indicating that the adhesion 
between the natural rubber and the Alloprene was weak. 

When the proportion of Alloprene was reduced to about 5% it was 
found that the dark spots were smaller and of lower contrast (fig. 9). The 
proportions are similar to those of the neoprene/natural rubber system of 
fig. 4, and it may well be that in both cases the proportion of the minor 
constituent was sufficiently small to render the polymers completely 
miscible. In such a case the Alloprene or neoprene molecules would be 
randomly kinked and interwoven with the molecules of the natural rubber 
as are the molecules in the bulk of a pure rubber (see Treloar 1958). The 
contrast would then be due to the greater local concentration of chlorine 
and would be low, and the edges of the darker regions would be expected 
to be diffuse. Experimentally, it was found that the contrast was very 
low and that the darker patches which often had diffuse edges could only 
be seen in regions where the specimen was very thin. The contrast could 
be increased by operating the electron microscope at 40 kv instead of the 
usual 80 kv but was still low. 

However, whether or not the polymers are miscible, there is one point 
of particular interest. The smaller dark patches are less than 100A in 
diameter and are of the order of size of the individual polymer molecules. 
It may well be that the smaller patches contain only one more or less tightly 
curled molecule. Information on this point was obtained by using a 
fractionated Alloprene, having a narrow range of molecular weights which ~ 
was prepared for the author by the British Rubber Producers’ Research 
Association. Films of purified natural rubber containing 5% of this 
material were prepared and examined and the size distribution of the dark 
spots measured. A histogram showing the size distribution of over 
1000 spots is shown in fig. 10. Spots greater than 180 A in diameter were 
not measured since they are outside the range of interest, nor were spots 
less than 404 in diameter since they become indistinguishable from the 
granular fine structure of the films. This fractionated Alloprene was 
stated to have a mean molecular weight of 70000 + 10000 (determined by 
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viscosity measurements). This would correspond to a sphere having the 
bulk density of Alloprene (1:6) of diameter 50+ 104. Thus the peak in 
the histogram at 80A probably corresponds to the diameter of a single 
molecule of the Alloprene, the slightly greater size being due either to it 
having a non-spherical envelope or to it being interwoven with the 
surrounding rubber molecules. The smaller spots presumably correspond 
to the tail of the molecular weight distribution. 


Fig. 10 
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SPOT DIAMETER (2) 


Histogram showing the distribution of spot size in natural rubber films 
containing about 5% Alloprene. 


It would be interesting to study the appearance and distribution of the 
spots in stretched films and an attempt to do this has been made 
Unfortunately the results were not reproducible. The film often stretched 
by thinning similar to that shown in fig. 3, but a few micrographs show 
evidence of the spots themselves. Figure 5 is an example. This film was 
not stretched deliberately but was deformed slightly by the stresses set 
up on drying and the spots appear elongated. 
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ABSTRACT 


Measurements have been made of the thermal conductivity of germanium, 
silicon, pure (n=3~x 10%) indium arsenide and impure (n=10'*) indium 
arsenide. In each case the electronic contribution xe to thermal conductivity 
has been calculated and the phonon contribution xz estimated by subtraction 
of xe from the experimental values. The thermal resistance of the lattice has 
not generally been found to be proportional to the absolute temperature as 
expected for three-phonon processes but the trend is towards a higher power 
relationship. The lattice resistance of heavily doped indium arsenide has been 
found to be greater than that of the pure material. 


§ 1. IyTRoDUCTION 


THERE have been few measurements of the thermal conductivity of semi- 
conductors above room temperature. The possible application of semi- 
conducting materials to the thermoelectric generation of electrical power 
makes thermal conductivity a factor of increasing importance. At the 
same time it is of interest theoretically since if one also has knowledge of 
the electrical properties of the semiconductor the lattice component of 
the thermal conductivity can be deduced near and above the Debye 
temperature. 


§ 2. MerHop or MEASUREMENT 


Earlier experiments on indium antimonide (Stuckes 1957) indicated the 
need for increased accuracy. Although the same comparative method 
(Stuckes and Chasmar 1956) has been used, certain modifications in thermo- 
metry and in thermal contacts have been found necessary to improve the 
accuracy to +5% up to 450°c. 


2.1. Thermometry 


Many precautions must be taken to avoid errors in thermocouple 
measurements due to work hardening and temperature ageing effects. 


ft Communicated by the Author. 
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Continuous 0-005in. diameter thermocouples of T-advance (nickel— 
chromium and copper-nickel alloys) have been found to give the least 
trouble amongst the more sensitive thermoelectric materials although it is 
necessary to have two sets of apparatus in use so that one can be kept for 
room temperature experiments only. 

Temperature measurements in a vacuum furnace were found to be 
inconsistent when the thermocouple junctions were electrically insulated 
from the copper blocks because the presence of exchange gas was essential 
to bring the thermocouples into temperature equilibrium. Rather than 
discard the differential method of measurement (which necessitates 
electrical msulation between thermocouples) experiments were carried out 
in an inert atmosphere. Calculation of the sideways heat loss from the 
column was made assuming a uniform temperature surround. This was 
achieved by placing a cylindrical copper shield between, the column and the 
cylindrical container; the temperature variation down the copper shield 
was found to be less than half a degree at 450°c. 


2.2. Thermal Contacts 


High accuracy in this method depends on there being virtually no 
contact resistance between any member of the stack and its thermocouple 
blocks. Since a good electrical contact is an indication of a good thermal 
contact, voltage probe measurements were taken along the column with a 
constant current flowing. A graph of voltage against distance should be a 
series of straight lines with sharp changes of slope at each contact—a step 
in the graph at the interface would indicate a high resistance at the contact. 
Precautions were necessary to eliminate thermoelectric e.m.f.’s when 
probing the surface of the semiconductor. 

Indium—mercury and indium-gallium alloys, which are both liquid at 
room temperature, were found to be unsatisfactory as intermediate 
contact layers above 200°c. The stack was therefore assembled by a 
soldering technique using pure indium after copper or nickel plating the 
semiconductors. Very good permanent contacts between the standards 
and the copper blocks were made by a silver-copper eutectic alloying 
process. 


2.3. The Comparative Standards 


In all experiments FH stainless steel has been used as a standard. The 
material was supplied by the Brown—Firth Research Laboratory and hada 
composition and heat treatment very similar to that investigated by 
Powell in 1936 and his values of thermal conductivity have been used 
throughout. 

As a subsidiary experiment the FH standards were compared with some 
Armco iron, which was part of a batch extensively measured by Powell 
(1934) and found to agree with Powell’s results to +3°% up to 425°c. 
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§ 3. EXPERIMENTAL RESULTS 
3.1. Germanium 

Only one sample of germanium was available at the time of these experi- 
ments and this was n-type and had a resistivity of about 50hmem at 
300°x. Unfortunately, electrical probe measurements indicated that it 
was very inhomogeneous, having a high resistance or barrier region near 
the centre; this may not matter since the electronic component of thermal 
conductivity is found to be very small. The results of two sets of thermal 
measurements are given in fig. 1. 


Fig. 1 
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Thermal conductivity of germanium and silicon. 


Four samples of germanium ranging in resistivity from 3o0hmcm to 
0-001 ohm cm have been measured by Abelés (1958). Although there was a 
wide spread in his experimental points, up to 700°K, the thermal con- 
ductivity of all four samples can be represented by the line shown in fig. 1. 
Above 700°K the very impure sample deviates from the others and has a 
lower thermal conductivity. The results of Ioffe and Ioffe (1954) and those 
of Kettel (1959) are also shown in this figure. 


3.2. Silicon 


Two samples of single crystal silicon have been measured, one ‘n’ type 
and one ‘p’ type, both having a room temperature resistivity of about 
3ohmem. Great difficulty was experienced in plating the silicon before 
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soldering. A nickel layer was deposited by an electroless plating method 
(Sullivan and Eigler 1957); from this good ohmic contacts were obtained 
so long as great care was taken. The results of thermal measurements are 
given in fig. 1. As the contacts deteriorated rapidly above 300°c no useful 
measurements have been made above that temperature. 

As far as it is known no other measurements on silicon at high tempera- 
ture have been published. The value of 0-30 cals(em°csec)— is in the 
middle of some rather divergent room temperature values. The American 
Metals Handbook (1948) quotes 0-20cals(cm°csec)-!, Ioffe and Ioffe 
(1954) give 0-26 cals (em°csec)-! whilst 0-35 cals (em°c sec)-! has recently 
been reported by Steele and Rosi (1958). 


3.3. Indium Arsenide 
One pure sample (x ~ 3 x 101°) and one sulphur-doped sample (n ~ 1019) 
of polycrystalline indium arsenide have been measured. Electrical 
conductivity and Hall coefficient up to 500°c (fig. 2) were measured on a 


specimen from each ingot cut immediately adjacent to that used in the 
thermal measurements. 


Fig. 2 
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Electrical conductivity of two samples of indium arsenide. 


From the ingot of pure material a piece 1-5cem long was initially used 
(much larger than the usual samples). After one set of oer ee nee 
ments this specimen was cut into two further pieces 1 cm an, 0- ie Oe 
and the experiment repeated on each. The results are given in fig. 3. 
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Only one sample of impure material was used but three sets of measure- 
ments, each on a freshly assembled stack, were made and the results shown 
in fig. 3. 

Bowers ef al. (1959) investigated two samples of indium arsenide, one 
similarly pure (n~5 x 10!6) and one rather less heavily doped (n~ 1038) 
than that measured in these experiments; their results are indicated in 
fig. 3. 
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§ 4. ELEcTRON ConTRIBUTION TO THERMAL CONDUCTIVITY 
4.1. Theoretical Discussion 


In the presence of a temperature gradient electrons and holes in a 
semiconductor will diffuse down the temperature gradient transporting 
kinetic energy from the hot to the cold end of the semiconductor. If on 
their way, an electron and a hole meet and recombine, they may give up 
their potential energy to the crystal lattice. This may be large compared 
to their kinetic energy and in some circumstances may be a very important 
mechanism of heat transfer; the terms ‘ambipolar diffusion’ or ‘ioniza- 
tion energy transport’ are sometimes applied to it. A further possible 
mechanism of heat conduction that has been postulated, but not so far 
definitely established, is exciton transport. An exciton is a related 
electron-hole pair with a weak binding energy; diffusing excitons could 


clearly carry thermal energy in a manner similar to free electron-hole 
pairs. 
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In addition to this electron contribution, heat is conducted through the 
semiconductor by vibrations of the crystal lattice (phonons) so that the 
thermal conductivity is given by 


Sas Kelectronic AP Klattice* 


Neglecting any contribution from excitons, the electronic thermal 
conductivity ke may be expressed in the form 


k\2 
re= (2) To hipaa ee eee eee) 


where o=on + op is the electrical conductivity at temperature T°K and is 
the sum of the contributions from electrons and holes. The factor A may 
be expressed as a sum of terms 


A = (Nee ae Ap ae Anp . ° . . . ° e (2) 


where the suffixes n, p and np refer to electrons, holes and electron-hole 
pairs respectively. The various terms involve Fermi—Dirac integrals 
whose form is related to the scattering mechanism in the semiconductor, 
while A,,, also includes a term involving the energy gap. 

This subject has been discussed by Price (1955), Ioffe (1956) and 
Madelung (1957). 

In the particular case where the relaxation times for electron and hole 
scattering are given by taxcH’™ and 7H’? respectively, assuming that 
classical statistics apply, that the scattering is elastic and that the states 
distribution in valence and conduction bands is parabolic, it can be shown. 
that 
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where bn= we 8, Op=rvpt3; 
Eg is the forbidden gap width at 7'°K and 
On =p ne, Op = PHpé- 


The value of v depends on the particular mechanism of scattering of 
electrons; for example, it can range from — 3 for acoustic scattering to + L 
for polar scattering above the optical mode temperature and to +3 for 
either impurity scattering or electron-hole scattering. 

The ambipolar contribution App will be most important when, on op 
that is for essentially intrinsic samples with comparable electron and hole 
mobilities. The value of H,/kT' is usually very much greater than 
vn +vp+5 and is dominant in this term. 
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For partial degeneracy the A’s may be evaluated in terms of Fermi-— 
Dirac integrals (see Madelung). In the limit of complete degeneracy the 
results are again simple. For example, if the electron system is completely 
degenerate while the hole system is non-degenerate (classical) then for all 
forms of electron scattering 


2 
= 55, Ap=(»p+ =, Me ee Le) 
Oy (ay 
A= non Je , Sn 3\° 6 
Anp= — 3 {ie + io tot a oe 


where €n is the Fermi energy. In general, degeneracy tends to increase 
the term in brackets but the effect is small until the electrons become 
highly degenerate. 

With sufficient knowledge of the electrical properties of a semiconductor 
xe can, therefore, be calculated. So long as xe is not too large compared 
to the measured value of thermal conductivity, the lattice component ‘x, 
can, then be deduced with reasonable accuracy. 


4.2. ke for Germanium and Silicon 


The germanium sample used in the experiments was fairly pure and by 
comparison, with some results of Dunlap (1957) would be intrinsic by 500°K 
and the silicon by 600°K. Electrical data from Morin and Maita (1954 a, b) 
have, therefore, been used to calculate xe for intrinsic germanium and 
silicon. The ratio of electron to hole mobility has been taken from the 
extrapolated high temperature curves. The temperature dependence of 
mobility in germanium is found to be 7'-'*® for electrons and 7'-*35 for 
holes and in silicon is 7'~** for electrons and 7'~*3 for holes. Although 
acoustic mode scattering would lead to a 7° dependence, nevertheless 
this has been assumed as a first approximation taking v= — }. 

The results are given, in fig. 4. _Ambipolar diffusion is dominant in both 
materials there being virtually no kinetic energy transport in silicon up to 
1000°x. By comparison with fig. 1 it can be seen that the electronic 
component in silicon is negligible and the experimental results represent 
the lattice conductivity. 

The predominance of ambipolar diffusion in germanium is confirmed by 
the results of Abelés. He finds that the thermal conductivity of his pure 
sample starts to increase above 700°K (the temperature at which it becomes 
intrinsic) while that of his impure sample continues to decrease until 
about 900°K (where it becomes intrinsic). At this temperature « for the 
pure material is 0-008 cals (em°c sec)~! greater than that of the impure; 
this is equal to the calculated value for ambipolar diffusion in intrinsic 
germanium at 900°K. 

Toffe (1956) has also calculated the ambipolar contribution in germanium ; 
he does not state what assumptions he has made but his results shown by 
the dotted line in fig. 4 seem somewhat high. 
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4.3. Ke for Indium Arsenide 


The calculation of «xe in indium arsenide is extremely complex, for not 
only has it been shown to be weakly degenerate but it also has a non- 
parabolic conduction band so that the simple formulae given in § 4.1 
do not apply. Moreover the scattering mechanism has not been, estab- 
lished. Ehrenreich (1957) has considered the problem of electron scatier- 
ing in the analogous compound InSb and finds that in the temperature 
range 200-500°K polar lattice scattering is predominant whereas at higher 
temperatures electron-hole scattering may become important. 


Fig. 4 
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In view of these facts no attempt has been made to use Ehrenreich’s 
analysis for the ambipolar transport term (denoted as A by him). The 
simple formulae (§ 4.1) have been. used, assuming classical statistics and 

olar scattering that is v=. : 
: In the calculation of Anp a value for the gap ralthes E, equal to 
0:5 ev-3-5 x 10-27 has been used. The value of Ey (#'z at 0°K) of 0-5 ev was 
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obtained from thermoelectric measurements by Antell ef al. (1956) while 
the temperature dependence was found from optical measurements by 
Oswald (1955). From a survey of the literature an average value for the 
ratio pp/un of 1/64 at 300°K has been obtained, this being only slightly 
dependent upon impurity concentration. Owing to a lack of reliable 
information on the temperature variation of both electron and hole 
mobility the value of 1/64 has been taken as constant over the whole 
temperature range. 

The kinetic energy and ambipolar contributions are shown in fig. 5 (@) 
while the total electronic component is given in fig. 5(b). It can be seen 
that in, indium arsenide the ambipolar diffusion is less than the kinetic 
energy transport. 

An attempt has been made to estimate the effect of variations in the 
assumptions used in, the above calculations. For example xe for two 
alternative scattering mechanisms is indicated by the dotted lines in fig. 
5(b). If, as seems likely by comparison, with InSb, electron hole scattering 
dominates at the higher temperatures then xe will be increased to the 
values shown, by the upper curve; should however, acoustic scattering 
predominate, «e would follow the lower curve. 

From Oswald’s optical measurements a value for the gap width 
EH,=0-48ev was obtained, surprisingly smaller than the value of 0-5ev 
from thermoelectric measurements. This would lead to a reduction in Ke 
of about 6%. 


Components of xe for Various Scattermg Mechanisms Allowing for 
Weak Degeneracy 


Kinetic energy transport 
Seattering (c/on) An 
mechanism 


Ambipolar transport dn 


It might be expected that the assumption of classical statistics in a 
material known to be weakly degenerate, would lead to an underestimate 
of xe. It has been shown by Antell ef al. that for a similar sample 
measured over the same temperature range, the value of »* (the Fermi 
energy in units of 7’) will at no temperature exceed 0, 1 or 2, according 
to whether acoustic, polar or electron scattering predominates. Both 
components of Ke (see table) have, therefore, been evaluated for these 
degrees of degeneracy using formulae given by Madelung (1957) in terms 
of Fermi—Dirac integrals, while the classical case (n*—+— co) has been 
evaluated from equations 3 and 4. 

From these figures, the maximum effect on the kinetic energy transport 
is an increase of about 8-5 which occurs for acoustic scattering. The 
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effect on the ambipolar contribution is an increase of about 12% at 
700°K where £)/k7’, Eg/kT and vp + § are all of the same order of magnitude 
—it will be less at lower temperatures and more at higher temperatures. 

The assumption of a constant mobility ratio has almost certainly led 
to an overestimate of the ambipolar contribution, since this ratio probably 
decreases as the temperature increases. If the mobility is a function of 
temperature of the form poc7’* then a (extrapolated from Folberth 
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et al. (1954) and Antell et al.) would be 2-3 for holes and between 1 and 1-5 
for electrons, leading to a value for zp/jin at 750°K of about 40% of that at 
300°K. Taking this into account and the effect of the lower value for the 
gap width, the lower estimate of the ambipolar diffusion would be about 
10% of the kinetic energy contribution for all scattering mechanisms. 
This contrasts with the statement of Bowers ef al. that up to 750°K 
the ambipolar term is less than 6% of the kinetic energy contribution. 
Unfortunately, they do not state what assumptions were made in 
their calculation. 

Estimation of the electronic transport in the impure sample is less com- 
plex. The ambipolar contribution will be negligible owing to the small 
number of intrinsic holes. At the same time the electrons are highly 
degenerate so that the electronic component is given, by 


2 /k\2 
«= 5 (2) ole 


The results are given in fig. 5 (0). 


§ 5. Lartice THERMAL CONDUCTIVITY 
5.1. Theoretical 


The theory of lattice conduction is somewhat complex and has not yet 
been, completely formulated: (cf. Berman (1953), Peierls (1955), Klemens 
(1956)). 

Above the Debye temperature © phonon-phonon U-processes pre- 
dominate and the following laws should be obeyed 


W,«T  three-phonon processes 
W,«T? four-phonon processes 


where W, (1/x,) is the thermal resistance. The first equation applies 
where the major phonon interaction is one involving three phonons, 
either two phonons combine to form a third or one phonon, splits up into 
two. The second equation holds where either one, two or three phonons 
interact to produce three, two or one phonon respectively. 

Well below ©, the resistance due to U-processes should decrease 
exponentially as the temperature is lowered until additional components of 
thermal resistance (such as phonon collisions with crystal imperfections, 
impurities, or crystal boundaries) swamp that due to interphonon 
collisions causing a rapid rise in W,. 

Even at high temperatures if there is a sufficiently large impurity 
concentration, there may be a contribution to the thermal resistance, due 
to phonon, collisions with the impurities, which is comparable to the 
resistance arising from interphonon collisions. According to Peierls, 
above the Debye temperature, the impurities will contribute an additional 
thermal resistance which is independent of temperature: this depends 
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on their number remaining constant and on their scattering cross section 
for a given phonon energy being constant, since (7’>@) the average 
phonon energy will be independent of temperature. 


5.2. Germanium and Silicon 


Since «e is found to be negligible in silicon up to 600°K the experimental 
results represent the lattice conductivity. kz; has been deduced for 
germanium by subtraction of xe from the measured values. The lattice 
resistance for each material is shown in fig. 6. There is a linear relation 
between W, and T up to 600°K for silicon (the full range of the experiment) 
but only up to 500°K for germanium. From 500°K to 700°K there is a 
marked increase in the temperature dependence. In neither case does the 
linear portion extrapolate through the origin. 


Fig. 6 
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Estimated lattice resistance of germanium and silicon. 


Low-temperature measurements by White and Woods (1956) and by 
Carruthers et al. (1957) have shown that the exponential dependence 
expected for U-process is not obtained in germanium or silicon, but the 
thermal resistance in each case is approximately proportional to the 
temperature. This has been explained (cf. Berman et al. (1956)) by the 
presence of isotopes which scatter the phonons. As the temperature is 
increased the thermal resistance due to U-processes will increase and 
ultimately become dominant. It is not clear at Whee stage this will 
happen, (@ for germanium is 360°K and for silicon is 658°K) and it is eee 
not surprising that over the linear part W,, is found to be Eee to 
the absolute temperature minus a constant. The results of Carruthers 
up to 300°K, which agree closely with those of White and Woods, are shown 
by the dotted lines in fig. 6 where they can be seen to be mutually consistent 
with the high-temperature results: the merging of the two curves is 
particularly good in the case of germanium. 
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The logarithmic plot (fig. 9) indicates that as the temperature is 
increased, the thermal resistance in germanium is tending to the square law 
dependence expected for four phonon processes. Since ke is only 6% of 
the measured value at 700°K it seems unlikely that a false assumption 
in calculating it would account for the departure of W, from linearity. 
Silicon does not show this effect ; in this case measurements were not made 
above the Debye temperature. 


5.3. Indium Arsenide 

The lattice conductivity for both samples of indium arsenide has been 
deduced using values of «e calculated in § 4.3 and the results are given in 
fig. 7; the deduced values of Bowers are indicated by the dotted lines. 
The most striking feature is the much lower lattice conductivity of the 
impure sample (n = 101°) compared to that of the pure sample (n= 3 x 101), 
contrasting with the fairly close agreement obtained by Bowers for both 
his samples (n =5 x 1016 and n= 1018). However, the two impure samples 
are not identical and there may be increased scattering of the phonons by 
impurities in the more heavily doped material. 


Fig. 7 
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Estimated lattice conductivity of indium arsenide. 


The full lines in fig. 8 represent the thermal resistance of the two samples 
assuming that the impure is completely degenerate and that for the pure 
material, classical statistics apply, polar scattering predominates and the 
mobility ratio is constant. The part shaded between, the dotted lines on 
either side of the lower full line represents the range of values for 
W,, (n=3 x 10") taking into account the uncertainty of electrical properties 
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discussed in, §4.3. The maximum values of W,, are obtained for electron- 
hole scattering, = 0-5 ev, a constant mobility ratio and weak degeneracy. 
The minimum values correspond to Hy =0-43 ev, a temperature dependent 
mobility ratio, acoustic scattering and no degeneracy. 

The lattice resistance of the impure material could be reduced if Ke 
has been, overestimated by the assumption of complete degeneracy to the 
extent shown by the shaded curve, the lower limit representing acoustic 
scattering and classical statistics. This limit is extremely unlikely since 
even if there is only partial degeneracy, polar, electron-hole or impurity 
scattering is the more probable scattering mechanism to which the upper 
limit is a better approximation. | 


Fig. 8 
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Estimated lattice resistance of indium arsenide. 


It would, therefore, appear that the thermal resistance of the impure 
sample is greater than that of the pure and the additional component 
increases as the temperature increases. 

Dr. J. M. Ziman and Mr. F. W. Sheard (private communication), 
however, have estimated that the sulphur-doped sample should have an 
additional component due to scattering of phonons by impurities of only 
about 0-5 cals (em°c sec) and this should be independent of temperature. 
They have suggested that the calculated value of xe may be too large 
(and hence too large a value for W,) when there is scattering of electrons 
by optical modes instead of acoustic modes. 

It is interesting to note that again the thermal resistance does not 
show a linear relationship with temperature over the whole range and 
the extrapolation of the lower part of the curves in, neither case passes 
through the origin. Bowers’ results between 700°K and 1000°K also 
indicate a more rapid dependence on temperature. From the logarithmic 


P.M, G 
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plot (fig. 9) the lattice resistance of the heavily doped sample can, be 
seen to be approaching a square law dependence. 


§ 6. CONCLUSIONS 
Taking into account the accuracy obtained in measurements on stain- 
less steel samples of known thermal conductivity and the reproducibility 
of repeated measurements on, semiconductors the experimental accuracy 
of measurement using the techniques described here is estimated as 
+5% up to 450°c. 


Fig. 9 


100 


50 


20 


100 200 500 1000 


Lattice resistance against temperature for germanium, silicon 
and indium arsenide. 


From fig. 9 it can be seen, that in no case is the thermal resistance 
directly proportional to temperature over the whole range as expected 
for three-phonon processes, although germanium approximates to this 
at the lower temperatures. The trend in other materials is to a higher 
power relationship as found for indium antimonide (Stuckes 1957). 
Since a 7° dependence would be expected for four-phonon, processes it 
may be concluded that higher order phonon processes are operative in 
these materials. 
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An, observation of particular interest is that the estimated lattice 
resistance of the highly impure InAs specimen is greater than that for 
the pure specimen, the difference being more marked at the higher 
temperatures. As a simple theoretical picture would suggest that 
impurity scattering should yield an, additional thermal resistance inde- 
pendent of temperature, this result may be of both theoretical interest 
and also of practical importance in thermoelectric applications. 
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CORRESPONDENCE 


Hall Effect in Liquid and Solid Mercury 


By P. W. Kunpatyi and N. E. Cusack 


Department of Physics, Birkbeck College, University of London 
[Received December 21, 1959] 


Unvit about 1930 there was some doubt as to whether liquid metals 
showed Hall effects. In 1914, Fenninger failed to find Hall voltages in 
liquid Hg, but he also found a zero effect in solid Hg and the latter result 
has been contradicted in at least two subsequent experiments. In 1931, 
Kikoin and Fakidow showed that a liquid Na-K alloy possessed a Hall 
coefficient approximately equal to the free electron value of — 1/nec, 
where 7 is the number of conduction electrons per unit volume. Never- 
theless they too concluded that liquid Hg has zero Hall effect. Taken in 
conjunction with later results for solid Hg (e.g. those of Serduke and 
Fisher 1932), this would mean that the Hall coefficient, R, changes 
considerably on melting. .It seemed desirable to check this result in 
the course of developing apparatus for investigating the Hall coefficient 
of liquid metals generally. 

The measurement was—by conventional d.c. techniques much as 
described by Lindberg (1952), the specimens being confined between flat 
glass plates whose separation was measured by an optical interference 
method. <A Hall voltage, of the order 1 wv, was detected and was found to 
be proportional to the magnetic field for constant current through the 
specimen, and vice versa. The Hall coefficient was the same for three 
specimens of thicknesses 4:17, 2:03 and 5:-5110-%cm. These facts 
indicate that a true Hall voltage was being measured, not seriously 
influenced by any magnetohydrodynamic effects. The mean of a large 
number of values of R was —74:6x10> +4% e.m.u. and this was 
independent of temperature between room temperature and the M.P. 
(—38:9°C). 

The best way of describing electron states in liquid metals is not yet 
clear. Using the terminology normal in solid state theory, the conditions 
that R should have the free electron value are, according to Cooper and 
Raimes, (i) the Fermi surface should be spherical and (ii) the relaxation 
time should be isotropic. It seems plausible to assume that both 
these conditions should hold in a liquid devoid of long-range order. 
On this view the Hall coefficient of liquid mercury should be — 77-4 x 10-5 
e.m.u. assuming two valence electrons per atom. This agrees quite well 
with the observed value. The results of Schulz (1957) on the optical 
‘properties of liquid Hg also fit very well with the free electron theory but 
the optical properties are not perfectly well established because the 
measurements of Hodgson (1959) do not confirm those of Schulz. 
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The values of & for solid Hg varied somewhat from specimen to speci- 
men as did those of Serduke and Fisher (1932). It is possible that this is due 
to the different ways in which the mercury crystallized each time it was 
frozen, for the specimens may have contained fairly large oriented regions 
and the Hall coefficient is probably anisotropic. The values of R, which 
were independent of temperature between —39 and —71°c. ranged from 
—87 to —101x10% e.m.u.; those of Serduke and Fisher from —72 to 

-88x10-° e.m.u. There is, of course, no reason to believe that either 
of Cooper and Raimes’s conditions are fulfilled in the solid. It is interest- 
ing that R is negative in the solid Hg while it is positive in Zn and Cd. 
This is a possible point against Wallace’s (1955) theory of bonding in these 
metals which treats them all on the same footing and entails the positive 
R in Zn and Cd. 

The present results for liquid Hg disagree with those of Kikoin and 
Fakidow. If the value of R now reported is correct, then, other things 
being equal, the Hall voltage from Hg will be about 17° of that from 
the Na-K alloy used by Kikoin and Fakidow. Since the error reported 
in their value of R(Na—-K) is + 10%, it is possible that their method was 
not sensitive enough to detect reliably the effect with Hg. Nevertheless 
it would be desirable to repeat the measurement with an entirely different 
technique, as well as to extend it to other metals, and such work is in 
progress. 
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The Effect of Electropolishing on the Retention of Vacancies 
During Quenching in Al and Al Alloys} 


By E. J. Fretset, M. E. Frne and A. Kerry 


Department of Materials Science, 
Northwestern University, Evanston, Illinois, U.S.A. 


[Received December 28, 1959] 


RecentLy Smallman and Westmacott (1959) observed diffuse scattering 
of x-rays at small angles in quenched pure metal foils. They attributed 
the effect to double Bragg scattering caused by distortion due to dislo- 
cations formed by condensation, of vacancies during quenching. Doherty 
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(a) Scattering curves for an Al-6 at. % Ag alloy: I—Electropolished and 
etched ; I1—Electropolished only. (6) Scattering curves for pure 
aluminium : I—Electropolished and etched ; I1—Electropolished only. 
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and Davis (1959) noted that surface oxide on aluminium prevented 
formation of surface pits which were due to condensation of vacancies, 
Herein are reported small angle scattering studies of foils (25 to 50 
microns thick) of Al and an Al-6at.°% Ag alloy quenched from 540°c 
indicating that a surface film is effective in increasing the magnitude of 
the double Bragg scattering ; escape of vacancies to the surface is reduced. 

The small angle scattering was measured using a General Electric 
XRD-5 diffractometer using a balanced filter technique to isolate the 
scattering due to CuK, radiation. The readings obtained were corrected 
for parasitic scattering. The foils were all capsulated in evacuated 
pyrex ampules (10-*mm Hg pressure) and drop quenched in water from 
540°C. 

The figure gives the intensity of the small angle scattering for two 
conditions: I—Electropolished (in a bath of 817c¢m? orthophosphoric 
acid, 40cm? water, 156g chromic acid, and 143cm®* sulphuric acid) and 
etched in hydrofluoric-nitric acid mixture; and Il—Electropolished 
only. It is to be noted that the small angle scattering below 1° is 
considerably larger for the unetched samples. Sample foils of Al-6 at. % 
Ag with the as rolled surface gave curves essentially identical to fig. 
1(a)—I. In the case of the alloy the scattering peak between 1° and 2° 
is due to the presence of Guinier—Preston zones. 

The following interpretation of the data in fig. 1 is suggested. In thin 
foils the surface is a very important sink for quenched-in vacancies. 
During electropolishing a surface film is presumed to form on the specimen, 
reducing escape of vacancies. The additional quenched-in vacancies 
then result in an additional number of dislocation loops or spirals. Htching 
removes the film and the effect is not observed. 
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REVIEWS OF BOOKS 


Llectrophoresis. Theory, Methods and Applications. Edited by Mian Brer. 
(Academic Press Inc., 1959.) [Pp. 563.] 107s. 6d. 
ELecrropuorssis has long played a major role in the study of colloidal systems 
and this volume sets out to give an account of the various forms of electro- 
phoresis which have been employed and the systems which have been studied 
by this technique. The first three chapters are largely theoretical and are 
devoted to consideration of the distribution of electric potential at interfaces, 
acid-base: equilibria of proteins and the theory of the moving boundary method 
of electrophoresis. Other chapters deal with the methods and applications of 
moving boundary, paper, zone and microscopic electrophoresis. The use of 
electrophoresis as a technique for the separation of mixtures and as a means of 
studying interacting systems is also discussed. The material, in general, appears 
to be soundly presented and well illustrated. The text is relatively free of 
typographical errors, but nomenclature, unfortunately, is not kept quite con- 
sistent from chapter to chapter. The emphasis of the book is very much in the 
direction of biological chemistry, a predominance not suggested by the title, and 
it is to be regretted that a little more space has not been devoted to the appli- 
cations of electrophoresis to other large and important fields, such as, mineral 
dressing, hydrophobic sols, emulsions, and suspensions in non-aqueous media. 
This book does, however, make a useful and commendable contribution to the 
subject of electrophoresis. 
R. H. O. 


Semi-conductors. By R. A. Smita. (Cambridge University Press, 1959.) 
[Pp. 494+ xvii.] 65s. 


THE output of work on solid state physics is now so great that the task of keeping 
up-to-date is an extremely difficult one. It is true that there are various series 
of summarizing reviews, but many of the articles are written before the subjects 
they deal with have been properly consolidated, and the student is often forced 
to read. not only the reviews but the original papers as well. It is therefore a 
pleasure to find a book such as Dr. R. A. Smith’s Semi-conductors which deals 
with a well-defined section of solid state physics in a workmanlike fashion. 

The first part of the book gives an account of the fundamental principles of 
carrier concentrations in semi-conductors, of electron transport phenomena 
including the diffusion of electrons and positive holes, and of optical effects. 
The treatment is far from devoid of mathematics, but, since the various steps 
are spelled out in more detail than in most of the standard texts, it should be 
generally acceptable to physicists and engineers. The second part of the book 
deals with the methods of determining the characteristic properties of semi- 
conductors, followed by an account of the properties of elemental and compound 
semi-conductors and brief accounts of their more important applications. The 
chapters dealing with the properties of the individual semi-conductors are 
particularly valuable. 


A. H. W. 
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Fig. 2 
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0-5 
Sequence showing the annealing of loops in quenched aluminium. Stage 
temperature 187°C. The interval between exposures is ~1 min. 
Note the loop A which remains through the whole sequence. Note 
also the slip trace B left by a moving dislocation in (b) and (c). Mag. 
x 53 000. : 
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(g) 10008 (h) - 10) 
Sequence showing the annealing of an isolated loop near a dislocation A. 
Stage temperature 197°c. The numbers refer to the time in seconds 

at which exposures were made. The loop lies on an inclined plane 

and the more steeply inclined parts are less visible. Mag. x 82500. 
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Sequence showing the annealing of loops near irregular dislocations. Stage 


temperature and time not recorded. As the loops anneal out the irregular 
dislocations become smoother. This effect is thought to be due at 


a least partly to climb. Mag. x28 000. 
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As-quenched dislocation loops. 


Three types of dislocation loops: (A) single loop, (B) hexagon, 
(C) double loop. (Note crystallographic alignment of the loops.) 
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z Annealed 10 min at 125°c after quenching. 
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Fig. 5 


Annealed 10 min at 125°c after quenching. 
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Annealed 10 min at. 125°c after quenching. 
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Annealed 10 min at 175°c after quenching. 
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yCold rolled 5°% after quenching. 


Fig. 12 


Cold rolled 15% after quenching. 


R. VANDERVOORT and J. WASHBURN Phil. Mag. Ser. 8, Vol. 5, Pl. 13. 
Fig. 13 


Cold rolled 15% after quenching. 


Fig. 14 


Cold rolled 45°, after quenching. 
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Natural rubber film. Very little structure is visible. 


Fig. 2 


Film of natural rubber with 10° of carbon black (‘ MPC black ’) filler. 


Fig. 3 


2% of carbon black (* MPC black ’) filler, after 
being stretched. Numerous thin patches have developed in the film, 


Film of natural rubber with 
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Film of natural rubber containing about 5°% neoprene. The dark patches 
probably represent areas of neoprene, the smaller ones corresponding to 
individual molecules. 
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Fig. 5 


Film of neoprene containing about 5°% naturalrubber. This shows the expected 
appearance of light patches on a dark background. 
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Film of 50% natural rubber, 50° Alloprene. The two polymers have not 
mixed and the film is inhomogeneous. 


Film of 80% natural rubber, 20% Alloprene. The dark patches are thought to 
represent globules of Alloprene. 
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Film of 90% natural rubber, 10°% Alloprene, after being stretched. The film 
has torn at the dark patches and this suggests that the adhesion between 
the two polymers was weak. 


Film of 95% natural rubber, 5%, Alloprene. Small dark patches—p bably 
due to individual molecules of the Alloprene—are visible. 


